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Abstract 

The physics of ferromagnetic doped manganites, such as Lai-xCaxMnOs with 
X PS 0.2-0.4, is reviewed. The concept of double exchange is discussed within the 
general framework of itinerant electron magnetism. The new feature in this context is 
the coupling of electrons to local phonon modes. Emphasis is placed on the quantum 
nature of the phonons and the link with polaron physics. However it is stressed 
that the manganites fall in an intermediate coupling regime where standard small- 
polaron theory does not apply. The recently-developed many-body coherent potential 
approximation is able to deal with this situation and Green's recent application to the 
Holstein double-exchange model is described. Issues addressed include the nature of 
the basic electronic structure, the metal-insulator transition, a unification of colossal 
magnetoresistance, pressure effects and the isotope effect, pseudogaps in spectroscopy 
and the effect of electron-phonon coupling on spin waves. 
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1 Introduction 



Since the discovery of high temperature superconductivity about fifteen years ago it has 
become clear that a metallic state obtained by doping an antiferromagnetic insulating oxide 
is substantially different from that of ordinary metals. Despite enormous experimental 
and theoretical effort neither the superconductivity nor the normal state of this unusual 
type of metal is at all well understood. There is no generally accepted view about the 
mechanism of superconductivity in the cuprates. To gain a wider perspective it is clearly 
important to study other oxides which become metallic upon doping. In fact interest has 
recently revived in the manganites Ai-xA^MnOs (A=trivalent rare-earth ion, A'=divalent 
alkaline earth ion), which were first studied fifty years ago. An excellent history of this 
early work is given in the introduction of an extensive review of mixed- valence manganites 
by Coey et. al. These compounds are ferromagnetic in the metallic state, typically 
for X ~ 0.2-0.4, and the origin of this metallic ferromagnetism is much less problematic 
than that of superconductivity in the cuprates. However above the Curie temperature Tq 
many of the manganites become insulating, and this metal-insulator transition is unlike 
anything that occurs in ordinary ferromagnetic metals such as Fe, Co, Ni or Gd. There are 
competing theories of this behaviour, and of the detailed exchange mechanism responsible 
for the ferromagnetism. Also the importance of electron-phonon coupling, and the possible 
existence of small polarons, is hotly debated. These are some of the issues which are 
focussed upon in this review. 

Two of the most-studied manganites are Lai-xSr^MnOs (LSMO) and Lai_xCaxMn03 
(LCMO). The parent compound LaMnOs, with nominal valence La^"'"Mn^+03^, is an an- 
tiferromagnetic insulator. Mn^+ ions have four d electrons and in the cubic environment 
the t2g states, with d wavefunctions of the xy, yz, zx type, lie lower in energy than the 
eg (x^ — y"^, 3z^ — r^) ones. The spins of the four d electrons are aligned, by Hund's rule, 
so the Mn^+ ion adopts a t2geg configuration with spin 2. The effect of doping, which 
removes x electrons per Mn atom, depends on whether LaMnOs is a Mott-Hubbard or 
charge-transfer insulator, in the terminology of Zaanen et. al. ||2|. In the latter case holes 
are created in the oxygen 2p band and, since the observed ferromagnetic saturation mo- 
ment is less than 4/iB per Mn atom, their spins must be antiparallel to the S = 2 Mn ions. 
In the former case electrons are removed from a narrow eg band and the three t2g electrons 
may be regarded as forming a local spin with S = 3/2. This is the traditional view @ 
which is still most widely held. The eg band contains n = 1—x electrons per atom with the 
possibility of metallic behaviour. The electrons are strongly correlated, avoiding double 
occupation of a Mn site, so that for n = 1 (x = 0) the system is a Mott insulator. The 
itinerant eg electrons are completely spin-polarized and align the local spins via Zener's § 
double-exchange mechanism. The validity of this general picture is thoroughly discussed 
in sections || and ^j. 

The structure of manganites such as LaMnOs is based on the ideal cubic perovskite 
(CaTiOa) structure shown in figure |I|. However the oxygen octahedra surrounding each 
Mn atom undergo two types of distortion. One is a rotation and tilting due to the La^"*" ion 
being smaller than the 0^~ and the other is a tetragonal distortion due to the Jahn- Teller 
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(JT) effect. Since tfie octaliedra are connected tfie rotations and tilting alternate in sense 
so that the final orthorombic structure has a unit cell whose volume corresponds to four 
perovskite unit cells. The JT distortion is predominantly of the Q2 type (see figure 0) 
in which each octahedron is stretched along one of the axes in the basal plane. Clearly, 
the locally stretched axis alternates in direction by 90°. The JT distortion arises from the 
degenerate Cg configuration; it splits the eg level so as to lower the energy of the Mn 3z^— 
orbital, with local z axis along the long Mn-0 bonds. The JT distortion in LaMnOs is a 
large effect, with an octahedral aspect ratio of 1.12 compared with 1.004 in CaTiOs 0. 
It encourages orbital ordering in the basal plane with local — orbitals pointing in 
alternate directions. This has been observed @. The magnetic ordering of LaMnOs, as 
indicated in figure |l|, is A-type antiferromagnetic, with ferromagnetic basal planes arranged 
antiferromagnetically. In a doped system such as Lao. 7810. sMnOs, which is a ferromagnetic 
metal, the crystal structure is much more nearly cubic f^. The slight distortion seems to 
be predominantly due to octahedral rotation and there is no evidence of a static JT effect 
Q. However a dynamic JT effect may occur as the eg electron moves from site to site and 
this constitutes a JT polaron. An important parameter in all manganites is the Mn-O-Mn 
bond angle which is less than 180° in a distorted structure. This reduced angle leads to 
reduced Mn-Mn hopping and hence to a narrowing of the eg band. 



Figure 1: The ideal cubic perovskite structure, showing the AF type A order for LaMnOs. 
Also shown are the breathing mode (Qi), the basal-plane distortion mode (Q2), and the 
octahedral stretching mode (Qs) which are present in the actual structure, (from reference 



Phase diagrams in the x-T plane of LSMO and LCMO are shown in figure ||. For 
X < 0.5 the two diagrams appear to be quite similar with a maximum Curie temperature 
Tc at X ~ 0.33. However the maximum Tq in LCMO is considerably lower than in LSMO 
and the paramagnetic state above Tq is insulating in LCMO and a poor metal in LSMO. 
In LCMO charge-ordering sets in at x ~ 0.5 and in the ferromagnetic insulator regimes 




i 



4 



at x < 0.2, in both LSMO and LCMO, substitutional disorder must play an important 
role. Much of the theoretical work reviewed here is concentrated on x ~ 0.3, where 
neither charge order nor orbital order is expected, and it is assumed that the system 
is homogeneous. However, it has been argued by Moreo et. al. and Nagaev |^ that 
charge inhomogeneity is widespread in the manganites with a tendency for holes, produced 
by doping, to segregate in ferromagnetic clusters. Unfortunately most quantitative work 
along these lines has employed models in which long range Coulomb forces are neglected. 
This leads to unrealistic predictions of macroscopic phase separation 0. 
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Figure 2: Phase diagram for Lai-xSixMnOa (a) and Lai-xCaxMnOs (b). The various states 
are: paramagnetic insulating (PI), paramagnetic metal (PM), canted insulating (CI), fer- 
romagnetic insulating (FI), ferromagnetic metal (FM), antiferromagnetic insulating (AFI) 
and charge-ordered insulating (COI). Tq and Tn are Curie and Neel temperatures, respec- 
tively. 

The importance of the Mn-O-Mn bond angle was mentioned above. This can be 
varied in the system Ao.yAg gMnOa, with a fixed number of electrons n = 0.7 in the Cg 
conduction band, by varying either the average A-site ionic radius {ta) or by applied 
pressure. Rodriguez-Martinez and Attfield ||ll| have also stressed the role of the variance 
in ta. An increase of either (r^) or applied pressure tends to push the Mn-O-Mn bond 
angle closer to 180°, thereby increasing the Mn-Mn hopping parameter and consequently 
the width of the eg band. The sensitivity of the Curie temperature Tq to such changes is 
shown in figure |^. The effect on the transport properties is even more striking. The contrast 
between the temperature-dependence of the resistivity p{T) for LSMO and LCMO, with 
X ~ 0.3, is seen in figures || and ^ In LSMO the resistivity above Tq (~ 370 K) continues 
to rise, as in a poor metal, whereas in LCMO the resistivity peaks at Tc (~ 260 K) and 
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decreases as the temperature is raised above Tq. This indicates a metal-insulator transition 
at Tc. Furthermore in LSMO, above Tc, p ~ 4 mf2cm whereas the resistivity of LCMO 
peaks at about 40 mf2cm . A satisfactory theory of the manganites must be able to account 
for this huge difference in behaviour between two apparently very similar materials. The 
effect of applied magnetic fields of magnitude 0-5.5 T on p{T) in LCMO is shown in figure ^ 
The large change in resistance near Tc is termed 'colossal magnetoresistance' (CMR) and 
it is this phenomenon which has inspired much of the recent research on the manganites. 
However the large fields required make it unlikely that this intrinsic property will be used 
in a sophisticated device for sensing magnetic fields. Low-field magnetoresistance observed 
in polycrystalline La2/3Sri/3Mn03 is attributed to spin-polarized tunneling between 'half- 
metallic' grains ||12|, |l^]. The half-metallic property of manganites is discussed in the next 
section on electronic structure. 
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Figure 3: The phase diagram of Ao.rA'QgMnOs as a function of the bandwidth. Closed 
circles represent variations of the bandwidth due to internal pressure (variations of the 
average A-site ionic radius (r^)), and open circles indicate variations due to externally 
applied pressure, (from reference [|T^] ) 



2 Electronic structure 

Before discussing the origin of metallic ferromagnetism in the doped manganites it is nec- 
essary to understand their electronic structure. The most important question is whether 
doping LaMnOs introduces holes predominantly into a Mn d band or into an O p band. 
This depends of course on which occupied band lies higher in energy. The question must 
be answered by a combination of theory and experiment. The usual starting-point is a 
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Figure 4: Resistivity versus temperature for Lai-xSixMnOs crystals. Arrows indicate the 
Curie temperature Tq for the ferromagnetic phase transition, (from reference [pTS]] ) 
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Figure 5: Resistivity (solid lines) of Lao.esCao.ssMnOa, acquired in applied fields ranging 
from to 5.5 T, and the bulk magnetization in an in-plane applied field of 50 G (•), both 
function of temperature, (from reference ||16|| ) 
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band calculation based on the local spin-density approximation (LSDA). In this approx- 
imation all electrons of a given spin experience the same crystal potential and therefore 
all d orbitals of that spin, apart from a modest crystal field splitting between t2g and eg, 
have the same on-site energy. In order to obtain a d"^ configuration, with n ^ 0, b or 
10, it is inevitable that the d band must overlap the Fermi level. Consequently in LSDA 
calculations for LaMnOs, with a Mn d^ configuration, the d band lies above the filled O p 
band ||, |l8l. The JT lattice distortion of Q2 type is required to open a small gap in 



the eg band so that the system becomes an insulator with orbital ordering. Also, the JT 
distortion stabilizes the observed antiferromagnetic structure over the ferromagnetic one 
H. The calculated Mn magnetic moment is 3.4/iB compared with 3.7/iB found by neutron 



scattering ||18|. Thus the LSDA calculations are compatible with the observed magnetic 
and structural details in LaMnOs. 



However it is well-known that in some oxides, like NiO |T9|, the p band lies above 
the d band and the LSDA fails ||2^. This is the case of a charge-transfer insulator where 
excitations across the gap are O p — > Ni d. To describe such a system within band theory 
one must use extensions of the LSDA, such as LDA+U or the self-interaction correction 
(SIC), or a Hartree-Fock method. In these approaches different d orbitals may see different 
potentials. For example one ] spin eg orbital could lie well below the Fermi energy and 
the other | spin eg orbital could lie above it. Since a particular eg orbital must be picked 
out for occupation on a given site cubic symmetry is broken and orbital ordering occurs. 
The splitting between the Cg levels would be the Hubbard on-site Coulomb interaction 
U. The p band can then lie above the occupied eg band. The density of states (DOS) 
obtained by Su et. al. pT| in a Hartree-Fock calculation for LaMnOs is shown in figure 0. 



Here the occupied majority spin d states (t2geg) lie well below the oxygen band and the 
Mn moment is close to 4/iB- The DOS sketched by Satpathy et. al. ||^ to illustrate the 
results of a LDA+U calculation for LaMnOs is slightly less extreme. The majority spin d 
band lies just within the bottom of the oxygen p band. LaMnOs orders magnetically as 
an A-type antiferromagnet (AAF) with ferromagnetic ordering within the basal planes. In 
figure P it is seen that the projected On — p DOS, corresponding to oxygen in the basal 
planes (see figure [l|), has an exchange splitting of about 0.75 eV at the top of the band. If 
the system were ferromagnetically ordered one would expect a similar, or slightly greater, 
exchange splitting. Thus, according to this picture, the ferromagnetic metal obtained by 
doping LaMnOs sufficiently with Sr, for example, would contain S = 2 Mn ions aligned by 
itinerant p holes of opposite spin. 

Both of the very different theoretical pictures presented above can explain the magnetic 
and structural properties of LaMnOs- To distinguish between them we must examine 



some spectroscopic data. Barman et. al. |^ give a very clear picture of the situation 



in the related material LaNiOs. This is a non-magnetic metal with the Ni^+ ion in a 
nominal t2geg configuration, that is ^IgGg^ for each spin. Thus the band structure should 
be very similar to the majority spin band in ferromagnetic Lao.sSro.sMnOs. Indeed the 
calculated DOS in reference ||22|, shown in figure 0, resembles in detail majority spin bands 
found in LSDA calculations for LaMnOs 0|, La2/sCai/sMnOs Ijl^] and La2/sBai/sMnOs 
23|| . Comparison should be made with the total DOS in figure ||, due to Pickett and Singh 



8 



6 



Mnd 



On p 



0l p » 4 



40 
30 
20 
10 



Total 



-12 



-10 



-6 -4 
Energy (eV) 



-2 



Figure 6: The (projected) DOS of LaMnOa, with AAF ordering. Positive and negative 
DOS are for up- and down-spin states, respectively. Energies are relative to the top of the 
valence band. The projected Mn-d and On-p DOS are for Mn and O on an up-spin basal 
plane, (from reference |[2]| ) 
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[TtII . Their calculation is for ordered La2/3Cai/3Mn03 with planes of La and Ca in a regular 
sequence. The feature of the DOS marked A in figure |^(b) is an eg band and B corresponds 
to the t2g band. C and D are two features of the broad oxygen p band. In the inset of 
figure ^(a) a calculated x-ray photoemission (XP) spectrum, including the effects of matrix 
elements as well as resolution and lifetime broadenings, is compared with the experimental 
spectrum. The peak near the Fermi energy Ep arises from the d band features A and B. 
The shoulder on the left of this peak and the second peak arise from the p band features 
C and D respectively. Sarma et. al. [l^ showed that this good agreement between theory 
based on LSDA band calculations and the observed XP spectrum extends to LaCoOs, 
LaFeOs and LaMnOs. In the Fe and Mn compounds, both insulators, the calculated 
valence band had to be shifted by 2.0 and 1.3 eV, respectively, to enlarge the band-gap. 
The LDA invariably underestimates band-gaps in semi-conductors and insulators (maybe 
LDA + small U should be considered, as discussed at the end of section |3T| ). In the 
ultraviolet photoemission (UP) spectrum of LaNiOs shown in figure 0(a) features A and 
B appear much more weakly than C and D, in contrast to the XP spectrum, and this is 
consistent with the assigned d and p angular momentum character |^. The UP spectrum 
of LaMnOs, both doped and undoped, is quite similar to that of LaNiOs with clear p 
band features C and D ||2j, However features A and B are not visible and seem only 
to contribute to a tail in the spectrum which extends to Ep in metallic samples such as 
Lai-xSixMnOs |^. It seems even more difficult to photo-excite d electrons near Ep, with 
ultraviolet photon energies, in doped LaMnOs than in LaNiOs. Furthermore only a weak 
shoulder corresponding to features A and B is seen in the XP spectrum of Lao.rSro.sMnOs, 
unlike the large peak in LaNiOs p^. However Park et. al. pBl, in high resolution spectra. 



see clear Fermi edges at 80 K in Lao.erCao.ssMnOs and Lao.rPbo.sMnOs. A contributory 
factor to low density of states at Ep may be the effect of electron-phonon coupling, which 
tends to open a pseudo-gap in the one-electron spectrum as discussed in sections | and p73| . 
Surface preparation is clearly also very important since the probing depth in photoemission 
may be as small as 5 A p7|| . 

According to LSDA calculations, for example that shown in figure ||, a doped manganite 
like La2/3Cai/sMn03 is almost 'half-metallic'. This term applies to a ferromagnetic metal 
in which the Fermi level lies in a gap in the DOS for states of one spin. In figure H the 
Fermi level lies just above a gap in the minority spin DOS. According to this calculation 
states from just below Ep to about 1.4 eV below it should be 100% spin polarized. Spin 
polarization at Ep can be measured directly by Andreev reflection and values of up to 80% 



have been found in Lao.rSro.sMnOs at 4.2 K by Soulen et. al. ||28| and Osofsky et. al. P9[ . 
Large spin polarization in Lao.eySro.ssMnOs is also deduced from spin-polarized tunnel- 
ing measurements ||30| , [3]| . Spin-polarized photoemission can probe below Ep and Park 
et. al. ||32| flnd 100% polarization in thin films of La^ySro sMnOs down to 0.6 eV below 
Ep. This suggests a smaller gap in the minority spin DOS than predicted theoretically. 
In early work on spin-polarized photoemission from Lao.7Pbo.3MnOs Alvarado et. al. USSj 
attributed their failure to see more than 20% polarization to a broad band of unpolarized 
Pb 6s states extending almost up to Ep. 

A stringent test of a proposed electronic structure of a metal is whether it gives the 
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Figure 7: (a) He I UP spectra of LaNi03 and NdNiOs. The experimental (dotted line) and 
calculated (solid line) XP spectra of LaNiOa are shown in the inset, (b) The total DOS 
(solid line), Ni d (dot-dashed line), and O p (dashed line) partial DOS for LaNiOs and 
NdNiOs calculated within LMTO-ASA H. 



observed Fermi surface. In a disordered system like Lao.ySro.sMnOs the Fermi surface is 
not sharp as in a pure metal. Even at T = there is no actual discontinuity in occupation 
number between states inside and outside a perfectly-defined surface in fc-space. Standard 
methods of determining the Fermi surface, like the de Haas-van Alphen method, can 
generally not be applied, except for the case of weak scattering in dilute alloys. However 
if scattering is not too strong the transition between occupied and unoccupied states in 
fc-space is sufficiently rapid to show up in the electron momentum distribution as a fairly 
clear Fermi surface. The momentum distribution can be measured by Compton scattering 
or positron annihilation. Livesay et. al. (3^ have used positron annihilation to investigate 
the Fermi surface of Lao.zSro.sMnOs and compare it with a LSDA band calculation. They 
used a cubic perovskite structure and a virtual crystal approximation, which makes the 
Fermi surface sharp, to find the band structure shown in figure ^. The system is almost 
half-metallic and the two majority spin Fermi surface sheets are shown in figure |10|. These 
comprise hole cuboids at the R {±n/a, ±n/a, ±7i/a) points that touch an electron spheroid, 
centred at the F (0,0,0) point, along the (111) directions. Similar Fermi surface sheets 
are described by Pickett and Singh |^. The analysis of the positron annihilation data 
is not simple, but Livesay et. al. conclude that their results agree well with the theory 
and establish the existence of the cuboid hole sheets. The virtual crystal approach is 
reasonable since Pickett and Singh ||17|, ^] have shown that scattering due to random 
A-site occupation has little effect on the majority-spin bands around Ep. However the 
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disorder localizes states at the bottom of the minority-spin band near Ep. Even if the 
bottom of the band falls just below Ep the system will still behave as half-metallic in 
transport ||35|. 
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Figure 8: The total DOS and local DOS on each inequivalent 
La2/3Cai/3Mn03 with FM spin ordering ||17 
planes sandwiching that layer of Mn ions. 
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The standard method of investigating band structure experimentally, away from the 
Fermi level, is angle resolved photoemission (ARPES). Mcllroy et. al. ||16| have investigated 
LaxCai-xMnOs and LaxBai-xMnOa with x = 0.35 using this technique. As pointed out 
above, the eg states close to the Fermi energy are hardly visible in photoemission, so one is 
looking at the dense ensemble of bands between 2 and 8 eV below Ep (see figure ||). It is 
hardly to be expected that any clear dispersion curves will emerge and Mcllroy et. al. deduce 
that in the Ca-doped compound there are a lot of very flat bands. They conclude that all 
the occupied states are extremely localized, but this seems an improbable assessment of 
the broad p bands. Some dispersion is found in the Ba compound, possibly associated with 
the steeply dropping bands in the TX direction at around 4 eV below Ep. A discussion 
of ARPES measurements by Dessau et. al. [3^ on a bilayer manganite is deferred to 



section 8.3.1. 
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Figure 9: Spin-polarized band structure of Lao.ySio.sMnOs. The majority spin bands are 



shown as solid lines, and the minority as dashed lines |^ . 
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Figure 10: Two sheets of the Fermi surface of Lao.rSro.sMnOs. Hole cuboids at the R 



points, coined 'woolsacks'. Electron spheroid centred at the F point |^ 
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Clearly there is a lot of evidence that the electronic structure of the manganites is well 
described by LSDA calculations and that states near the Fermi level have predominantly 
Mn d character with eg symmetry. However hybridization with O 2p states is by no means 
negligible. Sarma et. al. ||18| find that at the valence-band top and conduction-band bottom 
in LaMnOs the d component is 58% and 73% respectively. Simple one-band or two-band 
models must be thought of as hopping between orbitals of eg symmetry centred on Mn 
but extending onto neighbouring O atoms. The simple nature of the bands crossing the 
Fermi level, and the simple Fermi surface, suggest that a suitably parameterized two-band 
model, as used by many authors, may be quite realistic. Sarma et. al. [l^ discuss the 
reasons for the success of LSDA in the LaMOs compounds with M=Mn, Fe, Co and Ni. 
In these perovskites the M-0 bond is shorter than in the rocksalt MO compounds, so that 
the M-0 hopping parameter is larger. Also the effective Coulomb interaction U on the 
M site is smaller due to better screening by the closer O ions and to larger effective d-d 
hopping. 

We now briefly discuss some work which is often quoted in support of the hypothesis 
that LaMnOs is an insulator with a charge-transfer (O p — > Mn d) gap. Instead of a band 
calculation for the crystal Saitoh et. al. ||2^ performed configuration-interaction calcula- 
tions for an isolated MnOe cluster. The parameters of the model are the d-d Coulomb 
interaction U, the p — > d charge-transfer energy A and p — > d hopping parameter. They 
varied the parameters to obtain the best fit to the observed Mn 2p core-level photoemis- 
sion spectrum, O Is x-ray absorption spectrum and valence-band photoemission. They 
conclude that for LaMnOs U = 7.8 eV and A = 4.5 eV so that the condition U > A 
for a charge-transfer insulator is satisfied. On the other hand Chainani et. al. in a 
similar, but slightly different, analysis find A = 5.0 eV, U = 4.0 eV which makes LaMnOs 
more Mott-Hubbard-like than charge-transfer-like. Even for spectroscopic properties the 
results of cluster calculations seem inconclusive. Energy loss spectroscopy (EELS) on 
Lai_xSrxMn03 in which electrons are excited from O Is states to empty p-like states, has 
been presented [^] as confirmation of the hypothesis that LaMnOs is a charge-transfer- 
type insulator. A peak at the Fermi level grows with increasing x and this is interpreted 
as proof that doping produces holes in the oxygen p band. However it seems that the sce- 
nario extensively discussed in this section, where the charge carriers have Mn eg character, 
although significantly hybridized with O 2p states, is not excluded. 

In the next section we discuss the exchange mechanism responsible for ferromagnetism 
in the doped manganites from the general viewpoint of itinerant electron magnetism. It is 
clearly related to the nature of the underlying electronic structure. 



3 Itinerant electron ferromagnetism and double exchange 
3.1 The double-exchange (DE) model 

All work on the electronic structure of the manganites is in agreement on one point. The 
three majority spin t2g orbitals are all occupied and the minority spin ones are unoccupied. 
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This suggests that t2g electrons play no part in transport and may be considered as local- 
ized. Consequently the system may be considered as consisting of local spins = 3/2 on 
each Mn site coupled to electrons in a predominantly eg conduction band by local exchange 
interactions. According to the alternative view of the electronic structure, discussed in sec- 
tion 0, one eg electron per Mn site is also localized and local spins S = 2 are coupled to 
holes in the oxygen p band. Although we mention this alternative view again later we shall 
primarily adapt the former view for which strong evidence has been given in section |l]. 

Metallic rare-earth materials can also be considered as systems of local moments coupled 
to electrons in a conduction band by local exchange interactions. We now discuss the 
distinction between these and the manganites. The Hamiltonian for such a system is 

H = J2 UAc,. -jY.S,-<T,-hY,{St + at) , (3.1) 

ija i i 

where c}^ creates an electron of spin a on lattice site i, Si is a local spin operator and 
CTj = (erf, erf, erf) is a conduction electron spin operator defined by 

c't = (^i + icr^ = , cri = (^'t - = cJiCiT ^ = \ (^^T " ^n) (3-2) 

with riia = cl^Cia- The three terms of equation (|3.1|) describe hopping of the conduction 
electrons, exchange coupling between local and itinerant spins and coupling to an external 
magnetic field. If the local exchange coupling arises from hybridization between the local- 
ized and itinerant electrons, as in anomalous rare earth systems exhibiting heavy fermion 
behaviour, the exchange parameter J is negative. The Hamiltonian ( p.l|) is then often 
called the Kondo lattice model in view of its connection with the Kondo impurity model 
which has a local spin on one site only |^]. When Hund's rule coupling is dominant J > 
and the system is sometimes called a ferromagnetic Kondo lattice. This is misleading since 
for J > there is no connection with the Kondo effect. 

For J > it is useful to distinguish two distinct physical regimes, depending on the 
magnitude of J compared with the width 2W of the conduction band. If J -C W, as in a 
normal rare earth metal, the exchange coupling can be treated as a perturbation which gives 
rise to the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between local moments. 
In most rare earth metals this interaction, which oscillates in space, leads to oscillatory or 
spiral configurations of the localized f electron moments. The uniform ferromagnet Gd is 
an exception. In this weak coupling regime the Hamiltonian (|3T|) is usually referred to as 
the s — f or s — d model. 

If J ^ W the exchange coupling can no longer be treated as a perturbation. A 
conduction electron can only hop onto a site with its spin parallel to the local moment 
at that site. Furthermore if the number of conduction electrons per atom n < 1 double 
occupation of a site is strongly suppressed. A single electron at a site, with its spin 
parallel to the local spin S, enjoys an exchange energy —JS/2 which is lost if a second 
electron hops on. The system is therefore a strongly correlated electron system, just like 
the Hubbard model in the regime of strong on-site Coulomb interaction U, and for n = 1 
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the system is a Mott insulator. In much of the theoretical work on the present model the 
local spins are treated as classical vectors, corresponding to S" oo. Since for J 3> the 
itinerant spin must always be parallel to the local spin on each site, the effective hopping 
integral for hopping between sites i and j becomes Uj cos {9ij/ 2), where % is the angle 
between the classical spins Si, Sj. The cosine factor arises from the scalar product of two 
spin 1/2 eigenstates with different axes of quantization. The resultant band narrowing in 
the paramagnetic state favours ferromagnetism in order to lower the kinetic energy. This 
mechanism for ferromagnetism was first introduced by Zener ||^ and developed by others 



40| , p| . Since it involves strong exchange coupling on two adjacent atoms it is known 
as double-exchange. Consequently the Hamiltonian (|3.1| ) in the strong-coupling regime 
J ^ is called the double-exchange (DE) model. We note that the p-hole picture is 
not in the DE regime since the antiferromagnetic coupling between the Mn local spins 
and the O p-holes is weak compared with the large oxygen band-width. This review is 
largely concerned with the DE model, with quantum and classical local spins, and with 
its extension to include coupling of the electrons to local phonons. We call this extended 
model the Holstein-DE model |^. We now proceed to place the DE model in the wider 
context of itinerant electron ferromagnetism. 



3.2 A wider view 

The macroscopic exchange energy in a ferromagnet takes the form 

A I ^(Vn,(r))MV (3.3) 

where ni, n2, are the direction cosines of the magnetization direction at position r. 
Here A is the Bloch wall stiffness constant. A long wavelength spin- wave is a macroscopic 
oscillation in which the magnetization precesses around the equilibrium z direction. The 
transverse components of magnetization vary in space and time as 

M+ = + iMy (X exp [i{q ■ r - cut)] . (3.4) 

The dispersion relation is 

hu = Dq^ (3.5) 

and the spin- wave stiffness constant D oc A/Mg where Ms is the saturation moment. This 
macroscopic picture is valid for any ferromagnet, at zero or finite T (< Tq) and with an 
ordered or disordered crystal structure. In exciting a spin- wave quantum (magnon), for 
example by neutron scattering, a single spin is fiipped down. For ferromagnets with a 
substantial moment per atom, of a Bohr magneton or two, long wavelength spin-waves 
make an important, and often dominant, contribution to the temperature dependence of 
the magnetization for T < Tc/2, through the well-known Bloch T^^^ law. In manganites 
with Tq > 350 K the ratio 5 = D/{k^TQa?), where a is the lattice constant, is close to 
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the value 0.286 for the spin 3/2 nearest neighbour simple cubic Heisenberg model. For 
manganites with lower Tq, S can be larger by up to a factor 2 and this is attributed to 



strong electron-phonon coupling [^]. However in general D provides a good measure of 
Tc- This is not the case, however, in very weak itinerant ferromagnets such as ZrZn2, with 
small saturation moment, where longitudinal spin fluctuations are a determining influence 
on Tc. 

A very useful feature of D is that it is in effect a ground state property. It is directly 
related to the Bloch wall stiffness constant A, which can be determined from the ground 
state energy of the system in a non-uniform transverse magnetic field. Thus in princi- 
ple D can be calculated exactly in spin density functional theory, even though the band 
structure used in the course of the calculations has no such rigorous validity. In practice 
one obtains good values for D in Fe and Ni, and their alloys, on the basis of LSDA band 
calculations ||4^, |4^. We are aware of only one calculation along these lines for any man- 
ganite material. Solovyev and Terakura ||46|| find D ^ 300 meVA^ for Lao.yBao.sMnOs 
which is comparable to experimental values of somewhat less than 200 meVA^ for similar 
manganites with x = 0.3. We discuss the general approach below since it gives insight 
into the exchange mechanism in ferromagnetic metals and makes clear the basis of the 
double-exchange model. 

Edwards and Muniz showed how D can be evaluated directly from a multi-band 
tight-binding parameterization of a LSDA band calculation for the ferromagnetic ground 
state. No input about the underlying electron-electron interaction is required. However at 
a more microscopic level one may consider the band structure to arise from a Hartree-Fock 
approximation (HFA) to a Hamiltonian of the form 

H=J2Y1 (^)^ fc. Vfc- + ^int = ^0 + H,,t , (3.6) 

k fifj,'a 

where c^^^^ creates an electron of spin a in a Bloch state of wave-vector k formed from 
orbitals /i. Here the first term Hq represents the electron 'kinetic energy', actually the 
kinetic energy plus a spin-independent local potential, and i^int is a local (intraatomic) 
interaction term. Off-diagonal elements V^^' describe hybridization between different or- 
bitals. In general i^mt contains screened Coulomb interactions and exchange interactions 
of the Hund-rule type between the various orbitals. It can also include an additional one- 
electron term representing spin-independent diagonal disorder and may be generalized to 
include interaction with local phonons, as we shall discuss later. The Hamiltonian may be 
considered at two levels; either with effective interactions designed to reproduce the LSDA 
band structure within the HFA, or as a true many-body Hamiltonian in which correlation 



effects must be treated explicitly, for example by dynamical mean field theory |^, In 



both approaches one may start with an exact expression for D 49 



D = (lim(;.g)-^[J-,^^J) - limx.(0,..)) . (3.7) 

Here = Xlfc^i '^|tfc+<ji'^MfcT ^ spin-flipping operator representing a transverse component 
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of spin density and is the spin current operator defined by 



hqJ- = [5-, H] = Y,Y. (^MM'(fc) - y,Ak + q)) Vfct • (3. 



N-^ and Ni are the total number of electrons of | and | spin and the spin current- 

spin current response function, is given by the Fourier transform of a retarded Green 
function in the form 

xM,oj) = j dt((J,(t),J^,))e— *. (3.9) 



The expression in brackets in equation ( p.7| ) is proportional to the Bloch wall stiffness 
constant A (oc D{N^ — Ni)) and the second term is a generalized superexchange term 
which always makes a negative contribution to D. The first term therefore yields an 
upper bound to D and corresponds to the variational ansatz S~ \F) for the state with 
a spin- wave of wave- vector q excited, where \F) is the ferromagnetic ground state. This 
zeroth approximation Dq is often the dominant term in itinerant electron ferromagnets with 
maximum spin alignment such as the doped manganites. It is therefore worth examining 
in detail and for a cubic crystal may be written in the form 

^0 = . \ . E E VV,,,(fc)(cU.Vfc.) , (3.10) 

^ ^ ^' k fifi'a 

where the expectation value is evaluated in the ferromagnetic ground state and differenti- 
ation in is with respect to fc. If a particular orbital fia is completely full or completely 
empty in the ground state it makes no contribution to Dq. This follows because under 
these circumstances (cj^fccrC^'fea) and (cj^/fco.'^/ifcff) both zero for fi fJ^' and the integral of 
V^V^^(fc) over the whole zone may be written as a surface integral over the zone boundary 
which clearly vanishes. Thus in the manganites the t2g orbitals do not contribute to the 
first term in the Bloch wall stiffness constant A, as long as their mixing into j spin states 
above Ep and into | spin states below Ep is negligible. Of course the t2g electrons make 
the largest contribution to the factor N-^ — in D, and without their exchange coupling 
to the eg electrons there would be no ferromagnetism. The considerations above provide 
a justification for treating the t2g electrons as local spins, despite the substantial width 
of the t2g bands, and for the model DE Hamiltonian of equation (|3.1| ). We shall see that 
recent calculations of D in this model are special cases of the general formalism presented 
here. 

Quijada et. al. |^ have considered a two-band generalization of equation ( ^.ID and 



calculate D in an approximation equivalent to the random phase approximation (RPA). 
Nearest-neighbour hopping integrals on a simple cubic lattice are parameterized appro- 
priately for Cg orbitals, assuming only one non-zero Slater- Koster parameter (ddu). For 
J oo the result is just Dq given by equation (|3.10|) , with iV| — iV| = (25* + n)N where 



n = 1 — X is the total number of electrons per atom in the conduction bands and is the 
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number of atoms in the crystal. For nearest-neighbour hopping V^V^^'(fc) = — a^V^^/(fc), 
where a is the lattice constant and the on-site orbital energy taken as zero. Thus 

where K = N^^{Ho) is the expectation value of the 'kinetic energy' per atom. The K 



defined in reference [SQI is one-third of that defined here and has the opposite sign. In 



general our i^' < and equation (|3.11|) emphasizes that exchange stiffness is driven by 
one-electron ('kinetic') energy in the DE model. 

It is important to notice that equations (|3.10|) and are still valid if terms corre- 

sponding to spin-independent diagonal disorder and local coupling to phonons are included 
in the Hamiltonian. This last point has been exploited to investigate the effect of electron- 



phonon coupling on D [4^, as discussed in section The result concerning disorder 
is important because it shows that the double-exchange mechanism for ferromagnetism 
does not depend on having mobile carriers. Localized states contribute to the one-electron 
energy just as well as extended ones, as long as the localization length is more than a few 
lattice constants. This presumably explains why Tq varies continuously through the FM- 
FI phase boundaries in figure ||. At low carrier density {x < 0.2) the occupied states are 
below a mobility edge, making the system an insulator, but double exchange still operates 
to give ferromagnetism. 

We have laid considerable stress on Dq as the most important, typically metallic, part 
of D. However in some systems the value ofD is subject to enormous cancellation between 
Do and the second 'superexchange' term. Thus in a very weak itinerant ferromagnet such 
as ZrZn2, with very small magnetization M, the bracket in equation (|3.7|) is proportional 
to M^, so that D oc M. In a rare-earth metal like Gd, or in an Anderson lattice model of 
a ferromagnetic heavy fermion system, it is a delicate matter in this formalism to derive 



the correct RKKY-like result |^. The situation would be similar in a calculation of D 



based on the electronic structure picture of tlg^l {S = 2) spins coupled via p holes. Zhao 
5^ criticizes the standard DE picture with S = 3/2 spins coupled by eg carriers. His 



criticism is based on a two-band estimate of D with finite J, which comes out very small 



or negative. His model is exactly that of Quijada et. al. ||5^ who estimate values of D 
in good agreement with experiment. Zhao does not estimate D in the alternative p-hole 
picture which he favours. It would be very surprising if the RKKY-like exchange of that 
picture were larger than that of the itinerant electron DE picture. It would be interesting 
to carry out proper calculations of D in the multi-band tight-binding formalism starting 
from a LSDA band structure (compatible with the DE picture) and from LSDA+U or 
Hartee-Fock band structure (the p-hole case) . It was pointed out earlier that Solovyev and 
Terakura ||4^ obtained a reasonable value for D for a x = 0.3 manganite by means of the 
LSDA. They calculated the energies of static spin spirals of different wave-vectors q and 
fitted them to a Heisenberg expression for the magnon energy with exchange extending 
over many neighbours (see e.g. |§3|)- ^'^^ small q this is equivalent to calculating the Bloch 
wall stiffness constant A which is exactly related to D. However there is no reason to place 
much reliance on the dispersion curve for larger q. Solovyev and Terakura had difficulty 
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reconciling their work with the observed doping dependence of D; experimentally in LSMO 



54j| D increases with x for 0.175 < x < 0.3, just as Tq does (see figure |). It is clear, 
from their analysis of a simplified two-band model, that their opposite conclusion is due 
to the double-exchange being largest at a; = since the eg band is then half-full. At a; = 
the system should be a Mott insulator. Probably the best calculational method would be 
LDSA+U with a U large enough to split the t spin eg band but not large enough to reverse 
the ordering of p and d states. A picture like this is sketched in figure 12 of reference 

25|| . This does not change the general picture of the electronic structure we presented in 



section ||. 



3.3 Spin-waves in the DE model 

The spin- wave stiffness constant D in the one-band DE model (equation (lOP ) is easily 
deduced from the general multi-band formulation of section |3.1| . The only role of the local 
spins is to provide an exchange splitting between t and | spin bands and to contribute 
25* to the local moment — N^. We can therefore use equation ( |3.7| ) and need an 
approximation to xj(0,co'). A reasonable approach is to make a local approximation in 
which the electron self-energy S is assumed to be a function of energy only, with no 
dependence on wave- vector k. This is formally exact in infinite dimensions |^ and is 
the basis of the coherent potential approximation (CPA) and dynamical mean field theory 
(DMFT). In this case the vertex correction in the spin current response function xj(0,u;) 
makes no contribution and xj can be evaluated using a simple bubble diagram which 
involves the product of two one-electron Green functions. This result is well-known from 
the more familiar charge current response function used to calculate electrical and optical 



conductivity [0, |4§|. Then equation ( p.7|) becomes |Q ^ ^ 
D 



(3.12) 



where efc = '^jtij exp[ik ■ {Ri — Rj)] is the band energy, Ri being the position of site i, 
/(e) is the Fermi function and the one-electron retarded Green function Gka for spin a is 
given by 

Gfc.(e) = [e-efc-S.(e)]"' . (3.13) 
The occupation number (^ifccr) is determined by 

M = -\j Gka{t)f{t)dt. (3.14) 

On substituting this into equation (|3.12|) , and applying Green's theorem to the sum over 
k in the first term, we find 

D = [QT:{N^~N{)]-'lm I Ae^\e)Y,[Gu^{e)-Gul{e)f\Veu?. (3.15) 
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This compact expression was previously obtained ||57| , pq] within the framework of a CPA- 
RPA calculation for ferromagnetic alloys. However it is clearly exact within any local 
approximation, for example DMFT, and can be used in the presence of coupling to local 



phonons. It shows how the cancellation between the two terms of equation (|3.12|) takes 



place as — 0. An equivalent result has been given by Lichtenstein and Katsnelson 

||53|| . A convenient form of this result is obtained by introducing the function M[E) ||59| 
defined by the integral 

M{E) = (^]/87^3) J \\/ek\dS (3.16) 
over the surface = E in k space. Here fl is the volume of the unit cell. Then at T = 

de / dEM{E)lm [GE^ie) - GEi{e)f (3.17) 

oo J 



where GEai^) is defined by equation ( |3.13| ) with replaced by E and rio- is the number 



of a spin electrons per atom in the band in the ferromagnetic ground state. Normally, we 
consider J large enough such that rii = 0, so that n-^ = n = 1 — x. 

If we make the Hartree-Fock approximation (HFA) to Gka we have S| = —JS/2, 
Sj^ = JS/2 and it is easy to show that equation ( |3.12D yields the RPA result of Wang |^ 
and Furukawa |^ for D. The second term involves a factor A^^ where A = JS is the 
exchange splitting in the conduction band. If a Hubbard on-site Coulomb term U ^ • nj|nj| 
is added to the Hamiltonian of equation ( p.l| ) the only effect on D, assuming the spins are 
already aligned completely, is that A is increased to JS + Un. If we then set S* = we 
recover the well-known RPA result for D in the Hubbard model (see e.g. [§^]). If J — > oo 
in the DE model the RPA result for D reduces to the first term, which is called Do in the 
previous section. If we go beyond the HFA for G^i, in the state of complete spin alignment, 
we find that even for J cxo (but finite S) there is always some | spin spectral weight 
near the Fermi level which may even cause the totally spin-aligned state to be unstable. 



at least for S = 1/2 0, ^]. Thus half-metallicity may be destroyed and, even if not, 
the Fermi level lies at the top of the half-metallic gap. This effect was known long ago in 
the Hubbard model and is due to processes in the self-energy S| in which a | spin electron 
flips its spin and appears as a | spin electron near Ep plus a magnon ^7\, |6^, |69| . 
The t spin electron and the magnon may bind to form low-lying | spin quasi-particle 
states. The magnon dispersion curve necessarily runs into a continuum involving these 
states for wave- vector q close to an | spin Fermi wave- vector ||69|. Kaplan et. al. ||7^ have 



carried out some numerical studies of this effect in one dimension but seem unaware of the 
earlier analytical work. When the self-energy Sj^ discussed above is introduced into G^i in 
equation ( |3.15| ), the effect on D is to add a negative term to Dq (the RPA value for J — > oo) 
which is attributed to magnon-electron scattering |^ ^ The magnon excites t spin 
electron-hole pairs. This has been investigated diagrammatically from the point of view of 
al/S expansion by Golosov ||7^ and the effect can reduce D in the two-dimensional DE 
model by a factor of 2. Golosov also discusses magnon damping due to magnon-electron 
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scattering; he finds that the inverse magnon hfetime is proportional to the sixth power of 
the wave- vector q, as expected for a three-dimensional itinerant electron ferromagnet with 
complete spin alignment in the ground state ||73| . 

Another approach to this problem is a variational one similar to the ones of references 
[6^, |6^. An obvious ansatz for the state with one spin- wave excited from the ferromagnetic 



state of complete spin alignment is 



(S, + J2 A'^pS.+p-A^^Pl) \F) (3.18) 
fcp 

where S~ is the total spin-fiipping operator for local and itinerant spins. It is difficult to 
choose the coefficients A^p optimally, but a convenient choice is 

Afcp = [Nin + 2S)]-\ep - 6p+,)(efc - + fkuX-k)'' ' (3-19) 

where 

^ = [N{n + 2S)]-' 5^(n«T)(ej+, - q) (3.20) 
I 

is the RPA spin- wave energy for J ^ oo. Then we find an upper bound on _D, for J — >■ oo, 
in the form 

D = D 1 V I^^pI' 

^ (Vep • Vy) (v+p-fc - ep' - ep + efc) ) 

Here the sums over p, p' are restricted to occupied states in the t spin band and the sum 



over k is restricted to unoccupied ones. This result agrees with that of Golosov ||72| to order 



The last term in equation ( p. 211) , involving Vep ■ Vep/ with p 7^ p' in general, does 



not fit into the exact form of D in the local approximation (equation ( p.l2| ). This is not 
unexpected because the | spin self-energy including electron-magnon scattering depends 
on wave- vector k, not only on energy. The last term in equation (|3.21|) corresponds to a 



vertex correction first discussed by Edwards [71 . 



4 Magnetism and transport in the DE model 

In the last section we saw that much is known about the low-lying excitations of the DE 
model. Recent work on the properties of the model over the whole temperature range 
often deals with the limit of classical local spins (5* 00). The most complete work along 
these lines is that of Furukawa using dynamical mean field theory (DMFT) which is 
the best local approximation. No analytical progress can be made within DMFT for the 
quantum-spin DE model and numerical calculations have not yet appeared. However there 
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is no doubt that quantum spins introduce effects wfiicli do not exist for classical spins. For 
example the low-lying | spin spectral weight discussed in the last section has total weight 
(1 — n)/{2S + 1) for large J (see section |4.2|) and therefore does not appear for — oo. In 
this section we shall therefore concentrate on another local approximation, the many-body 
dynamical CPA p4|, |75| in which quantum spins can be treated analytically and which 



agrees with DMFT in the classical limit. Like much current research on the DE model (see 
section the method has its roots in work done on the Hubbard model in the sixties, 
in this case by Hubbard himself ||76|. Since such work is no longer common knowledge we 
start with an introduction to CPA in the Hubbard model. 



4.1 CPA for the Hubbard model 

To introduce the many-body CPA we consider the Hubbard model, which is a simpler 
model for strongly correlated electrons than the DE model with quantum spins. The 
Hamiltonian for this model is 



and Hubbard set out to find an approximation to the one-electron retarded Green 
function Gfc(j(e) which is exact in the atomic limit tij = 0. He used the equation of 
motion method and the idea of the alloy analogy described below. It turns out that 
Hubbard's approach, without the minor "resonance broadening correction", is equivalent 



to the CPA which was developed later [^]. The CPA derivation of Hubbard's result is 
much simpler than the original equation of motion method. However we had to resort to 
an extension of the original method to derive the many-body CPA for the DE model, with 
and without phonons, in the general form needed to discuss magnetic properties. In this 
paper we restrict the derivation to the paramagnetic state in zero magnetic field, although 
we summarize some more general results. 

The alloy analogy consists in considering the |-spin electrons, say, to move in the 
potential due to static |-spin electrons, frozen in a random configuration which must be 
averaged over. Thus a one-electron Hamiltonian for |-spin is obtained from by taking 
the last term to be U ^■ni^{nii) where (riii) takes the value 1 with probability and 
with probability 1 — nj. Here is the number of cr-spin electrons per atom. It is important 
to note that the alloy analogy is quite distinct from the Hartree-Fock approximation in 
which {uifj) = Tifj for all i. In the alloy analogy a cr-spin electron moves in a random 
potential given by U on sites and on 1 — n^^ sites. In the CPA the random potential 
is replaced by a uniform, but energy-dependent, effective potential ^^(e) for an "effective 
medium". This effective potential, in general complex, is called the coherent potential 
and is in fact the electron self-energy. The procedure for determining Eg. is to insist on a 
zero average t-matrix for scattering by a central atom, with potential U or 0, set in the 
effective medium. Equivalently the average of the site-diagonal element ^^(e) of the Green 
function, for each type of central atom, is put equal to the site-diagonal element of the 
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Green function for the effective medium. Thus 



G„ = n. 



+ (1 -^a) 



G„ 



1 + KGo 



and 



Go (e-S.(e)) 



(4.2) 



(4.3) 



where the bare band Green function is given by 



de 



e — e 



(4.4) 



Here = X]j^ji6xp[ifc ■ (i?j — Rj)] is the band energy, where Ri is the position of site 
2, iVo(e) is the corresponding density of states per atom and N is the number of lattice 
sites. Equations (|4.2|) and ( ^1.31) are to be solved self-consistently for ^^(e) and hence for 
the local Green function G^{e). Equation ( |4.2| ) may be written as 



Gfj 



+ 



(4.5) 



This may be compared with the exact Green function for the atomic limit {tij = 0) which 



is given by ^6 



G 



1-ns- 



where in this retarded Green function e has a small positive imaginary part as usual. Hence 

G.(e) = Gf + . (4.7) 

Clearly this CPA equation is exact in the atomic limit, when A^o(e) = ^{^) and it follows 
from equation ([4.3|) that T,^ + G~^ = e. Solution of the CPA equation becomes simple if 
the density of states iVo(e) is taken to be of the elliptic form 



2\l/2 



where 2W is the bandwidth. Then from equation (0| 



(4.8) 



"e_(e2-iy2^i/2j _ ^^ ^^ 

Introducing this expression for in equation (|4.3| ), and solving for e — ^^(e), we find 

J:^ + G-^ = e-W^G^/A. (4.10) 



Go(e) 
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Hence equations ([4 .61 ) and ( |4.7] ) give an algebraic equation for Ga. 

Solving this type of equation for G^, Hubbard |^ calculated the density of states 
Nfj{e) = — 7r~^Im G^{e), considering particularly the paramagnetic state = = n/2 
and concentrating on the half-filled band case n = 1. He showed that for U/W greater 
than a critical value, equal to 1 in the present approximation, a gap opens in A^(e) at the 
Fermi level so that the system becomes an insulator as envisaged by Mott. For U ^ W the 
density of states consists of two peaks centred on e = and e = U, these being broadened 
versions of the 5-functions at these energies in the atomic limit. Furthermore, for general 
band-filling n, the spectral weights in the two peaks are the same as in the atomic limit. 

The CPA for the Hubbard model has some serious defects. There are no self-consistent 
solutions with magnetic order. Furthermore in the paramagnetic metallic state, for n < 1 
or for n = 1 with U/W less than the critical value, the system is never a Fermi liquid. 
There is never a sharp Fermi surface at T = with a Migdal discontinuity in the Bloch 
state occupation number, as pointed out by Edwards and Hewson |^]. This is due to the 
absence of states with infinite lifetime at the Fermi level, since within the alloy analogy 
all states are scattered by disorder. A modification of the CPA to remedy this defect, 



retaining the analytic simplicity of the method, had some limited success ||7^, |80 



However the most satisfactory approach is DMFT which involves numerical solution of an 



associated self-consistent impurity problem |47, 82|. DMFT may be regarded as the best 



local approximation, in which the self-energy is a function of energy only, and is exact in 
infinite dimension. 

The many-body CPA is considerably more satisfactory for the DE model than it is 
for the Hubbard model, as discussed in the next section. There is one limit, the case 
of classical spins {S = oo), in which the CPA is identical to DMFT. This is because 
classical spins are static and an alloy analogy of frozen disordered spins is completely 
justified. DMFT for the DE model has only been implemented fully for classical spins 
||83| , Isij and the many-body CPA discussed in the next section provides an approximate 
analytic extension of DMFT to quantum spins. The system orders ferromagnetically below 
a Curie temperature Tc, as it should, and the disordered spin state above Tq should be 
well described. However, the accuracy of the ground state at T = for finite S is unclear. 
The saturated state with all itinerant and local spins completely aligned, which is the 
ground state for S" = oo (we are always considering large J in the DE model), is never 



a self-consistent CPA solution for finite S \74 . This is due to low-lying | spin states 



of the type discussed in section |3.3| . In the treatment described there involving magnon 
emission the local approximation corresponds to neglecting the magnon energy, in which 
case I spin quasi-particle states always appear below the Fermi level and the saturated 
state is unstable. In their approximate treatments of DMFT Meyer et. al. |^ also find | 
spin states below the Fermi level for any finite S. Possibly this is a general result of the 
local approximation. The true parameter range of stability of the saturated ground state 
is unknown. It has been shown rigorously that for J = cx) it is unstable for S = 1/2 and 
0.12 < n < 0.45, with a simple cubic nearest-neighbour tight-binding band ||63|. Wurth 
and Miiller-Hartmann ||65| find no instability with respect to a single spin-fiip for any n 
when J = oo and S = 3/2. If the true ground state is not saturated it seems unlikely to be 
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a uniform (spatially homogeneous) ferromagnet, with partially ordered local and itinerant 
spins, as in the uniform CPA ground state for finite S. Such a state would probably not be 
a Fermi liquid, just as in CPA, unless the electrons making up the spin S became partially 
delocalised with spectral weight at the Fermi level. Numerical DMFT results for quantum 
spins are urgently required to throw light on the matter. 



4.2 The many-body CPA for the DE model 

Edwards et. al. fT^ developed the many-body CPA for the DE model using an extension 
of Hubbard's equation of motion method. Hubbard's "scattering correction" becomes more 
complicated owing to the form of the interaction term in the DE model whereby electrons 
can flip their spin via exchange of angular momentum with the local spins. This dynamical 
effect couples the equations for and G| and was flrst treated by Kubo |^ in a one- 
electron dynamical CPA. The main feature of the many-body CPA is that we recover 
Kubo's one-electron CPA as n ^ and the correct atomic limit for general band-filling n 
as tij ^0. A second paper ||75| showed the equivalence to DMFT in the limit S ^ oo, 
J ^ oo. The full equation of motion derivation of the many-body CPA is required to obtain 
general results in the presence of a magnetic field and/or magnetic order ||7^. However it 
turns out that in the zero-field paramagnetic state we can deduce the CPA equation from 
the atomic limit Green function G^^ and equation (|4.7|) , just as in the Hubbard model. 
We shall therefore not repeat the full derivation in this paper although we shall discuss 



results on magnetic properties in section 



The atomic limit Green function, G^^, say, is easily obtained by the equation of motion 
method using the Hamiltonian (|3.1| ) with tij = |7^. The result for zero field {h = 0) is 
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S^){l-n^)-S-a+) 



e + J{S + l)/2 



^ {(S + l + Sn{l-n^) + S-a+) ^ as 



e- J(5 + l)/2 
^1- S')n^ + S-a+) 



JS/2 



e - JS/2 



(4.11) 



and for /i 7^ one merely has to replace e by e + h/2. The angle brackets (. . .) represent 
thermal averages and all operators within them correspond to the same site i, this suffix 
thus being omitted. This expression, with four poles, is considerably more complicated 
than the two-pole Hubbard model expression of equation (|4.6|) . The poles at e = ±JS/2, 
±J{S + l)/2 correspond to energies to add or remove an electron from the atom, that 
is to transitions between singly-occupied states and either unoccupied or doubly-occupied 
states. The singly-occupied states have total spin 5 + | or S* — | with energies — JS/2 and 
J{S -|- l)/2 respectively; the unoccupied and doubly-occupied states have zero energy. In 
a state of complete spin alignment, with all local and itinerant spins t, Gf^ has a single 
pole at e = —JS/2 and Gf^ has three poles at e = ±JS/2 and J{S + l)/2. The weight in 
the low-lying | spin level at — JS/2 is (1 — n)/{2S + 1). 

In the zero-field paramagnetic case it turns out that the CPA equation for G{e) with 
the redundant suffix a omitted, is given by equation (|4.7| ) as in the Hubbard model. Thus, 
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taking the band to have the eUiptic form (|4 

G{e) = G^^ [ 



the CPA equation for G is 
- W^G/A) 



(4.12) 



with G^^ given by 



2S 



nS/2- {S-a) ^ 5(1 



n/2) - (5 ■ cr) 



_e + J{S + l)/2 e- J{S + l)/2 
{S + 1)(1 - n/2) + (5 ■ cr) n{S + l)/2 + {S ■ a) 



+ 



(4.13) 



e+JS/2 ' e-JS/2 

The spin symmetry of the paramagnetic state has been used to simphfy the expectations 
in the previous form of G^^, equation (|4.11j ). It is easy to show that {S ■ cr) nS/2 
as J — > cxD and {S ■ cr) will be very near this limit as long as JS > 2W. We make this 
approximation in calculating G(e), and hence the density of states, N{e) = — 7r~^Im G{e) 
from equations (|4.12|) and (|4.13|) . The results are shown in figure for S = 3/2 and 
J = AW for various n. Clearly, from equation ( |4.13|) , the approximation to {S ■ cr) has the 
effect of removing the weak band centred on e = —J{S + l)/2 but it does not affect the 
total weight or the distribution of weight between the two lower and two upper bands. It 
may be seen that as n increases from the band near e = J{S + l)/2 is reduced in weight 
and a new band appears near JS/2, until at n = 1 no weight remains in the band near 
J{S + l)/2. The weight in the band near —JS/2 is (S* + 1 — n/2)/{2S + 1) per spin so 
if JS is sufficiently large to separate the bands {JS > 2W) this band will just be filled 
at n = 1 producing a Mott insulator as expected. This redistribution of weight between 
bands as they fill with electrons is characteristic of the many-body CPA and was missing 
from Kubo's one-electron CPA |^ which was restricted to n = 0. 

In the strong-coupling limit J — > oo, which is taken with a shift of energy origin 
e — s> e — JS/2, equation ( ^.13|) simplifies to 

G^^{e) = t-\S + 1 - n/2)/{2S + 1) . (4.14) 



Equation (|4.12|) then becomes a quadratic equation for G with solution 

2 



Ve 



(4.15) 



where a'^ = {S + 1 — n/2) / {2S + 1) and D = aW. By comparing with equations ([4.8|) and 
( W^) we see that the density of states is a single elliptical band of weight and bandwidth 
2aW. As S" — i> oo the band-narrowing factor a 1/V2 = 0.707, which is close to the 
classical result of 2/3, obtained by averaging cos (6^/2) over the solid angle. 



In the classical spin limit S 
equation becomes 

1 
2 



oo we rescale J, replacing it by J/S, and the CPA 



G 



Here we have used the general equation 



+ 



J/2 S + G-i - J/2_ 



(4.16) 



(|4.7| ), valid for arbitrary band-shape, rather than 
equation (|4.12| ). Equation (|4.16|) is precisely the equation obtained by Furukawa [^] within 
DMFT. 
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Figure 11: The density of states in the paramagnetic state of the double-exchange model 
for S = 3/2, J = AW, and n = 0, 0.25, 0.75 and 1. The Fermi level is always within the 
lowest band until this is just filled for n = 1. Energy units of W are used, (from reference 

0) 
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4.3 Resistivity in the paramagnetic state of the DE model 

The Kubo formula for the conductivity o involves the two-particle current-current response 
function. However in the local approximation of CPA or DMFT there is no vertex correc- 



tion ||5^, ^ and a may be expressed in terms of the one-particle spectral function 

^fc(e) = -TT-^Im G'fc(e) = -Tr-^Im [e - - S(e)]-^ . (4.17) 

In the paramagnetic state G is T-independent if we assume {S ■ cr) = nS/2, and a de- 
pends on temperature only weakly through the Fermi function. If we neglect this thermal 
smearing around the Fermi energy /i we may calculate at T = but consider the results 
to apply to the actual paramagnetic state at T > Tq. We find ||7| 



^-^.E^llM,)? (4.18) 

k 

where Vk = Ve^ is the electron velocity and is the volume of the unit cell. Just as in 
section we can introduce the function M{E), defined by equation ( p.l6| ), and write 



^=^J dEM{E)\AM\\ (4.19) 

where Ae{^) is defined by the right-hand expression in equation ( |4.17| ) with e = ^, tk = E. 
For a simple cubic tight-binding band = —2t'^p cos k^a, with /3 summed over x, y, z, 
V^efc = — a^efc. Then it is straight-forward to show that 

/E 
EN^{E)dE (4.20) 
-oo 

where Nc{E) is the density of states for the simple cubic band. If Nc{E) is replaced 
by a suitably-scaled Gaussian Ng{E) = (3/7r)^/^l^~^ exp [—3{E/W)'^], corresponding to a 
hypercubic lattice in infinite dimensions, M{E) becomes Mg = {a^W"^ /6)Ng{E). However 
if A^c(-E') is replaced by the elliptic density of states No{E) it becomes Mq{E) = [a^CW"^ — 
E^)/3]Nq{E). We refer to previous treatments of the M{E) factor in a when we discuss 
the optical conductivity in section p. 3. 2 . 



Since it is convenient to use A^o(e) to calculate the Green function and self-energy, as 
in the previous section, it is reasonable to evaluate a using equation (|4.19|) with M{E) = 



Mq{E). Edwards et. al. ||7J| took the strong-coupling limit J oo for simplicity and 
the results for the resistivity p = are plotted against band-filling n for various S in 
figure |12|. The lattice constant was taken as a = 5 A which is comparable with the Mn-Mn 



spacing in manganites. For J = oo the band-width W is the only energy-scale and, since 
the integral in equation (|4.19| ) is dimensionless when a factor is taken out, p does not 



depend on W. In fact in the DE regime JS > 2W the resistivity is almost independent 



of both J and W. It is seen in figure |12| that p diverges correctly at = 0, owing to 



the absence of carriers, and at n = 1 where the system becomes a Mott insulator. Well 
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away from these insulating limits p does not depend strongly on S, so that quantum spin 
effects are not very important. Furthermore p ^ 1 mficm over a wide range of band- 
filling, which is much smaller than observed in some manganites above Tc, for example 
in LCMO (see figure ||). This agrees with the conclusion of Millis et. al. |^ that the 
DE model, with electrons scattered purely by disordered local spins, cannot describe the 
physics of the manganites completely. Early work by Furukawa |^ seemed to point to 
another conclusion, although the DMFT is equivalent to the CPA approach. Edwards 
et. al. ||74| showed that the confusion arose from Furukawa's use of a convenient, but rather 
unreasonable, Lorentzian density of states. This makes the calculation of Ae{p) very 
simple and Furukawa effectively took M{E) in equation ( [4.19|) to be Lorentzian also. In 
later work Furukawa [^] used the elliptic band. Results of such calculations for the limit 
J — s> oo are shown in figure O and it is remarkable that p is at least an order of magnitude 



larger than one finds in figure 12 for the more reasonable elliptic band. 




0.4 0.6 

Filling (n) 



Figure 12: The zero field paramagnetic state resistivity p = versus band-filling n 
for the double-exchange model from reference Here J = oo, a = 5 A, and S = 

1/2, 1, 3/2, 2, 5/2 and oo, p increasing with S. The elliptical density of states and formula 
(|T|) with M{E) = Mo{E) are used (from reference 



4.4 Magnetism in the DE model 

As discussed at the beginning of section |4^ , the full equation of motion approach to many- 
body CPA is required to determine magnetic properties such as spin susceptibility x ^^nd 
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Curie temperature Tq. In reference [|7^] this involved a hierarchy of Green functions satis- 
fying 45* + 1 coupled algebraic equations for local spin S; only the S = 1/2 case was briefly 
discussed. In reference ||75|| a major simpliflcation was achieved by introducing generating 
Green functions which generate all the required Green functions by differentiation with 
respect to a parameter. The coupled equations are then replaced by a single flrst-order 
linear differential equation, the parameter being the independent variable, whose analytic 
solution yields the required CPA equations for the Green functions. The classical limit 
S = oo can then be taken and for J = oo the equations coincide with those of DMFT, 
which are only obtainable in the classical limit. The many-body CPA is therefore an an- 
alytic approximation to DMFT for arbitrary quantum spin 5* which becomes exact for 
S = oo. The many-body CPA also coincides with Kubo's |^ one-electron CPA in the 
limit n ^ where that is valid. 

To determine the magnetic properties one problem remains; the CPA and DMFT equa- 
tions contain one set of correlation functions ((S"^)™) which cannot be obtained directly 
from the Green functions. There is an indirect procedure for determining these correlation 
functions within CPA but it proves to be unsatisfactory, never yielding ferromagnetic solu- 
tions. However, for S* = oo, DMFT provides a way to calculate the probability distribution 
function P{S^), and hence ((S"^)™), and Green and Edwards used an empirical exten- 
sion of this for finite S. This extension guarantees that the spin susceptibility exhibits 
the correct Curie laws for band occupations n = and n = 1. Thus for n = we have 
a Curie law over the whole temperature range, corresponding to N independent spins S. 
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For n = 1, with J = oo, we have independent spins S + ^. For < n < 1 we find a finite 
Curie temperature Tc and some results |^ are shown in figure |l^. In figure |14(a)| Tq is 
plotted as a function of n for various S with J = oo, using the elliptic band. Clearly for 
finite S ferromagnetism is more stable for n > 0.5 than for n < 0.5, in agreement with the 
findings of Brunton and Edwards ||63|. For S = oo the result agrees closely with that of 
Furukawa In figure |14(b)| we see the effect on Tc for 5* = 1/2 of changing the bare 
band-shape from elliptic to simple cubic tight-binding. A dip in Tc occurs around n = 0.3 
which is the region where the ground state of the simple cubic DE model with S' = 1/2 is 
rigorously not one of complete spin alignment Ip3 . 




(a) (b) 

Figure 14: The Curie temperature k-oTc/W of the double-exchange model versus band- 
filling n for various calculated with J = oo using the elliptic band (a); the effect on 
Tc for = 1/2 of changing the elliptic band to the density of states for a simple cubic 
tight-binding band with nearest neighbour hopping (b). (from reference ||75| ) 



5 The metal-insulator transition 

In section ID we pointed out the wide range of behaviour in manganites. LSMO with x ~ 0.3 
has a metal-poor metal transition at Tc, whereas in LCMO the resistivity p(T) decreases 
with rising temperature above Tc (see figure 1). Moreover the resistivity near Tc is an 
order of magnitude larger in LCMO than in LSMO. This behaviour of LCMO, and many 
other manganites, is characterized as a metal-insulator (MI) transition. Furukawa has 
pointed out that the contrasting behaviour of LSMO is quite well described by the DE 



model . Furukawa's p(T) curve (figure 9 of reference ||6T|) is similar to that of figure |2T 
in this paper, which corresponds to weak electron-phonon coupling i.e. essentially to the 
pure DE model. A central problem in the theory of the manganites is to explain the 
origin of the MI transition in systems like LCMO and to explain why it does not occur in 
LSMO. Clearly the pure DE model is insufficient, as first emphasized by Millis et. al. |^ . 



However a common element of many theories of the MI transition is the narrowing of 
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the band in the paramagnetic state above Tq due to the DE mechanism. This narrowing 
of the band encourages any tendency to localize the electrons and produce insulating 
behaviour. However theories differ greatly as to what actually drives the localization. 
Millis et. al. |^, |8^ and Roder et. al. propose that the driving force is strong electron- 



phonon coupling associated with the dynamic Jahn-Teller effect. On the other hand Varma 
||90|| and Miiller-Hartmann and Dagotto |^ proposed that localization could occur due to 
a combination of magnetic disorder and non-magnetic A-site disorder. Nagaev's theory 



1^, |9j is also based on A-site disorder while Furukawa ||6T| argues in favour of phase 
separation models of LCMO with, for example, ferromagnetic and charge-ordered regions. 
We begin with a discussion of models based on disorder. 

5.1 The role of disorder 



Sheng et. al. ||9^ pursued a line similar to Varma's |£0|. They studied localization due to 
both non-magnetic randomness and the off-diagonal disorder associated with the classical 
spin DE model. As mentioned in section pTT| , the effective hopping integral is tij cos{9ij/2) 



where the angle 9ij between neighbouring spins varies randomly in the paramagnetic state. 
Sheng et. al. showed using scaling theory that this off-diagonal disorder alone can only local- 
ize a small fraction of electron states close to the band edges but fails to cause localization 
of the electron states at the Fermi level for the ferromagnetic regime with x = 0.2-0.5. 
Li et. al. ||9^ came to a similar conclusion using a transfer-matrix method. This is an 
important result because local approximations such as DMFT and CPA are not able to 
detect localization. Both groups further showed that to obtain a MI transition between the 
ferromagnetic and paramagnetic state, by means of the mobility edge moving through the 
Fermi level on introducing the off-diagonal disorder of the paramagnetic state, one needs 
large diagonal disorder with random site energies in the range {—W, W). Here 2W is the 
bandwidth in the ferromagnetic state and ~ 1 eV typically. The random energies at 
Mn sites in Lai-xCaxMnOs, say, arise from varying local environments of La^"*" and Ca^"*" 
ions. Neglecting metallic screening, the change in potential at a Mn site due to changing 
the charge on one of the neighbouring A sites by e is about 0.4 eV, assuming a local di- 
electric constant of 10. Pickett and Singh ||35| deduce from a band calculation on ordered 



La2/3Cai/3Mn03 that in the disordered case Mn site energies have a Gaussian-like distri- 
bution with full width at half maximum equal to 0.6 eV. This is much less than 2W so the 



condition for a MI transition is not fulfilled. Smolyaninova et. al. |95| searched for scal- 
ing behaviour of their resistivity measurements on LacerCacssMnOs and Ndo.rSro.sMnOs 
thin films; they concluded their results were incompatible with an Anderson localization 
transition. 

Nagaev's |^, ^ theory is rather different. He assumes the charge carriers are holes in the 
O p band so that DE does not apply. His idea relies on the fact that a charge disturbance in 
a metallic ferromagnet produces a magnetic response extending over a magnetic correlation 
length ^. In a state of complete spin alignment at T = 0, ^ is essentially the lattice spacing 
but ^ diverges at Tc- Nagaev therefore proposes that A-site disorder leads to extensive 
static variation in the exchange potential which could localize carriers, particularly for 
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T ^ Tc- An applied magnetic field reduces ^ and will lower the resistance as required 
for the CMR effect. Unfortunately, no quantitative work seems to have been done on the 
model which makes it difficult to asses. 

Diagonal disorder in the DE model with classical or Ising local spins has also been 
treated using CPA ||9^ or DMFT ||9^, |98|. Within these local approximations localization 
does not occur and the only possibility of a MI transition is for a gap to open in the 
density of states at the Fermi level. A random site potential is considered, within a one- 
band model, which takes the value A/2 on a fraction x of the sites, and — A/2 on the others. 
The idea is that for a suitable A a gap will open as the band narrows in the paramagnetic 
state. These models are faced with several problems. A more Gaussian-like distribution of 
site energies, as appropriate on Mn sites due to A-site disorder, would not produce a gap 
and, even with the binary distribution, one requires A ^ W, the half band-width, which 
is unrealistically large. Another difficulty is the position of the gap relative to the Fermi 
level. The band which should be split by disorder is the lowest one in the density of states 



shown in figure [Tl|. In the case of classical or Ising local spins the total weight in this band 
is one state per site (including both spins) and, when split by disorder, the weight below 
the gap is 1 — x. The gap is therefore positioned at the Fermi level as required for n = 1 — x. 



However for quantum spins S it is clear from section [4.2| that the corresponding weight 
is (1 — x){2S + 2 — n)/{2S + 1); thus the electron density n would have to be assigned 
arbitrarily a strange dependence OTvx,n = 5(1 — a;)/(5 — x) for S = 3/2, to place the Fermi 
level in the gap. 

The conclusion of this section is that disorder alone cannot account for the MI tran- 
sition in manganites with x ~ 0.3. The effect of A-site disorder is sufficiently weak for a 
virtual crystal approximation to be a reasonable starting-point for considering the effects 
of electron-phonon coupling. Of course if this coupling is so strong that narrow polaron 
bands appear then disorder will become important 0]. We discuss this in section |^. 

5.2 Electron-phonon coupling 

For sufficiently strong coupling to local phonon modes an electron produces a strong distor- 
tion of the lattice in its immediate neighbourhood. This distortion moves with the electron 
and the whole structure is called a small polaron. At low temperatures the small polaron 
occupies a narrow band of coherent Bloch-like states which narrows even further as the 
temperature rises. At high temperatures, typically larger than that corresponding to half 
the phonon energy, the polaron moves diffusively with an activation energy (see e.g. ||99| ). 
A condition for small polarons to exist is Ep > Wz'^^"^ [|100| , where is the polaron 



binding energy, W is the half-width of the bare electron band and z is the coordination 



number of the lattice. The idea proposed by Millis et. al. |g3 is that as the band narrows on 
passing into the paramagnetic state, by the DE effect, this criterion just becomes satisfied. 
Thus in the low-temperature ferromagnetic state, where the criterion is not satisfied, the 
system is metallic and electron-phonon interaction merely produces a small enhancement of 
the quasi-particle mass; the quasi-particles in this Fermi-liquid are sometimes called large 
polarons. As Tq is approached small polarons are formed and their motion is governed 
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by an activation energy related to Ep. This picture of the MI transition was developed 
in different ways by Millis et. al. ^ and Roder et. al. HI. Millis et. al. used DMFT 



||88| and Roder et. al. ^ ^ 

and treated the phonons classically, whereas Roder et. al. used the quantum-mechanical 
approach of small-polaron theory. This approach is discussed further in section |^. 

The electrons may couple to different local phonon modes, which may have the Qi, Q2, 
Q3 symmetries shown in figure |l]. The breathing mode Qi is the simplest but the Q2 JT 
mode is usually regarded as the most important and was emphasized by Millis et. al. |^ . 
To describe the JT coupling properly one should consider a two-band model with the correct 
Cg symmetry. This was done by Millis et. al. ||8^ and Zang et. al. [|101|| . Roder et. al. |^ 



and Zang et. al. ||101|| treat DE in a simple mean-field way, by means of a band-narrowing, 
but avoid double occupation of sites by effectively using spinless electrons. Millis et. al. |^ 



treat the Hund's rule coupling J but in a two-band model large J no longer produces a 
Mott insulator for n = 1, as it does in the one-band model. It is therefore necessary to 
introduce on-site Coulomb interaction [[102| , [103|| . Millis et. al. |^ do not do this so that 
for n = 1, which should correspond to the undoped insulator, they have a metal. Also, 
for weak electron-phonon interaction, they find Tq is largest for this value of n. This 
contrasts strongly with figure |14| where, in the one-band DE model, Tc = at n = and 
1 and has a maximum in between. Held and Vollhardt ||103|| , using DMFT, obtain very 



similar results to those of figure 14 for = cx) in a two-band model with strong on-site 



Coulomb interaction included. However the values for Tc are about twice those in figure |l4 



Nevertheless, if on-site Coulomb interaction is neglected for simplicity, it is sensible to use 
a one-band model to describe the manganites. In fact Green |^ developed the theory of 
such a model (the Holstein-DE model) in which electrons couple to quantum spins and 
quantum phonons. He used the many-body CPA method ||74|, |7^ which was discussed 
in section ^j. In the limit of classical spins and phonons this method is equivalent to the 
DMFT of Millis et. al. [8^. Green's work therefore bridges the gap between this classical 
phonon treatment and the polaron theories. We therefore describe it in some detail in the 
next section. Comparison between theory and experiment indicates that electron-phonon 
coupling in the manganites is just in the intermediate regime where small-polarons are on 
the verge of forming. At present the many-body CPA is the only theoretical method which 
can deal satisfactory with this cross-over regime and in section || we show how it can be 
used to tie together experimental data using many different techniques. 



6 The many-body CPA for the Holstein-DE model 



This section is based on Green's |^ recent study of the Holstein-DE model in which the 
electrons of the DE model couple to local phonons as in the Holstein treatment of small 
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polarons [104 |105|| . The Hamiltonian is 



I] a 



(6.1) 



The first three terms constitute the DE Hamiltonian of equation (^Tj) while the first, fourth 
and fifth terms form the Holstein model. Einstein phonons on site z, with energy uj and 
creation operator 6|, couple to the electron occupation number rii = ^^riia with coupling 
strength g. The electron-phonon coupling is of the form —g'^^UiXi, where Xi is the 
displacement of a shell of atoms surrounding site i, and in application to the manganites it 
may be regarded as an effective Jahn-Teller coupling. Previous studies of this model have 
either concentrated on coherent polaron bands, like Roder et. al. |^, or have treated the 



phonons classically |^ so that there are no polaron bands at all. The many-body CPA 
approach is able to encompass both aspects and to describe the crossover from quantum 
polarons to the classical picture as temperature and/or model parameters are varied. The 
relationship to previous theoretical work and to experimental studies of the manganites is 
discussed fully in section However we mention briefiy below some related work on the 
pure Holstein model, without coupling to local spins. 

Sumi considered the Holstein model with one electron in the band, first treating the 
phonons classically ||106|| and later quantum mechanically [|107]| . The classical case, with 
frozen displacements Xj, corresponds to a multicomponent alloy for which CPA is the best 
local approximation. In his dynamical CPA treatment of quantum phonons, Sumi ||107 



treated the one-site dynamics correctly and his work is completely equivalent to the more 
recent DMFT treatment of Ciuchi et. al. ||1U8|| . As a general rule dynamical CPA and 
DMFT are the same for one-electron problems. DMFT is the correct extension of CPA 
to the many-body problem of finite electron density but for the Holstein model, as for 
the DE model, it cannot be carried through analytically in the quantum case. Numerical 
work ||109|| applying DMFT to the Holstein model has been aimed mostly at understanding 
superconducting transition temperatures and charge-density-wave instabilities rather than 
the polaron physics with which we are mainly concerned. An unfortunate feature of the 
Holstein model for spin 1/2 electrons is that in a quantum treatment the true ground 
state for strong electron-phonon coupling consists of unphysical singlet bipolarons with 
two electrons bound on the same site. This problem does not occur in the one-band 
Holstein-DE model since strong coupling J to local spins prevents double occupation of 
sites, as pointed out earlier. It is also bypassed if the phonons are treated classically, as 
in the work of Millis et. al. [|110|| on the Holstein model. The Holstein model is more 
complicated than the DE model and it turns out that Green's many-body CPA no longer 
reduces to the correct one-electron dynamical CPA/DMFT ||107] , p.08|| as band-filling n ^ 0. 
Although correct in the atomic limit tij = 0, the theory is clearly cruder for the Holstein 
and Holstein-DE model than for the pure DE model. 

We start by deriving the Green function for the Holstein-DE model in the atomic limit. 
The Hamiltonian Hal in this limit is given by equation (|6.1D with the first term omitted 
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and with site indices and summation suppressed. We remove the electron-phonon couphng 
by the standard canonical transformation |^ H = e'^i^ALe"^ where s = ~{g/uj)n{b'^ — b). 
Under this transformation b ^b + {g/uj)n and the Hamiltonian separates into a fermionic 
and bosonic component: 

H = H{ + H^ (6.2) 



H{ = -JS-(T-h{S' + a')- {g'^/iu)n\ H^ = iub^b. (6.3) 

The transformation corresponds to a displacement of the equilibrium position of the phonon 
harmonic oscillator in the presence of an electron and the downward energy shift g^/ u; is a 
polaron binding energy which we write as Xu, where A = g'^ I<jJ^ ■ If two electrons occupy the 
site (n = 2), which will not occur for large J, the energy shift becomes ^g^ joJ^ corresponding 
to an on-site bipolaron. Writing out explicitly the thermal average in the definition of the 
one-particle retarded Green function we have 

Gf{t) = -ie{t){[c^{t),ci]^) 

Tr|e-/5^A. [c.(t),4] ) 

Tr{e-^^AL} ^ ' 

and the canonical transformation introduced above can be carried out within the traces, 
using the property of cyclic invariance. Thus Hal ^ ^"^cl and c„{t) becomes 

e^^*Xc,e-^^* (6.5) 

where X = exp [(7(6^ — 6)/c<j] . Using equation (|6.2|), we can write the traces in equa- 
tion ( |6.4| ) as products of fermionic and bosonic traces. Hence we find 

(t) = -ie{t) {{c^{t)ci)^F{t) + {clcAt)),F*{t)} (6.6) 

where F{t) = (X(t)X^)b and the thermal averages ( )f, ( )b correspond to the systems 
with Hamiltonians iff and i^b respectively. It may be shown that 

F{t) = e-^(26+i) [5(5 + i)]i/2 [^(^ ^ ^^/2)] } (6.7) 

where b = b{uj) = (e'^" — l) ^ is the Bose function with f3 = (/cbT")"^. The last factor is of 
the form exp (zcos0) which generates the modified Bessel functions lr{z): 

00 

exp(zcos(/.) = ^ I,(;z)e^"'^. (6.8) 

r=— 00 

To evaluate the fermionic averages we consider for simplicity the limit J — » 00 in zero field 
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Figure 15: One-electron spectra of the Holstein-DE model in the atomic limit at zero and 
very high temperature. They consist of delta-functions, with energy spacing u;, whose 
strength is indicated by the envelope curves. The plots are for the paramagnetic state with 
S = J = oo, h = 0, n = 0.5, u/W = 0.05 and g/W = 0.18, where W is a unit of energy 
later to be identified with the half- width of the electron band in the full Hamiltonian P2[ . 
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{h = 0). Then the last term in Hf may be written —{g^/u)n, since n = or 1 only, and 
this may be absorbed into the chemical potential which is finally determined to give the 
correct number of electrons n per atom. Thus H{ is just the DE Hamiltonian in the atomic 
limit and the sum of the two fermionic averages corresponds to the function G^^{t) whose 
Fourier transform is given by equation ( ^■14| ). It is easy to see that the first and second 
thermal averages in equation (|6.6|) take constant values (1 — n){S + 1)/{2S + 1) and n/2 
respectively. Hence, from equations (|6.6|) - (|6.8|) , we obtain the Fourier transform of G^^, 



with J ^ oo and h = 0, in the form 

I,{2A [b{u;){b{u;) + 1)]'/'} (2^ + l)f e^^^-/^ + + l)(l _ n)e--/^-/2 



{2S + l)exp{X[2b{u) + l]} e + ru 

(6.9) 



The density of states — Tr^^Im G^^{e) is shown in figure for the classical spin limit 
5 oo at quarter-filling n = 0.5. It consists of delta- function peaks separated in energy 
by uj and the envelope curves show the weight distribution at low and high temperature. 
W is an energy unit which, when we go beyond the atomic limit, will be the half-width 
of the itinerant electron band, as usual. The values adopted for the parameters uj /W and 
g/W relate to the manganites, as discussed in section ||. The symmetry of the spectrum 
about zero energy is due to the choice of filling n = 0.5; in general at T = the lower and 
upper 'bands' have weights n and 1 — n respectively. By counting weights it may be seen 
that for any n the chemical potential lies in the peak at e = 0, which has very small weight 
e~^/2 per spin. The shape of the envelope function at T = 0, with two maxima and very 
small values at the centre of the pseudogap between them, may be understood physically 
as follows. The delta-function at e = ru {r > 0) corresponds to an excitation from the 
ground state, with no electron and the undisplaced oscillator in its ground state, to a state 
with one electron and the displaced oscillator in its r^^ excited state. The strength of the 
delta-function depends on the square of the overlap integral between the displaced and 
undisplaced oscillator wave functions. Clearly this is very small for r = and goes through 
a maximum with increasing r as the normalized displaced wave function spreads out. At 
T = it is easily seen from equation (|6.9| ), using lr{z) ~ (z/2)l''l/|r|! for small that the 
weight of the delta-functions at e = ±ru is proportional to A'/r! . Hence the maxima in 
the envelope curve occur at e ~ ±Xuj, which is the polaron binding energy. 

We now turn to the Holstein-DE model with finite band-width. As for the DE model 
it is necessary to use the full equation of motion method to derive the many-body CPA 
in the presence of a magnetic field and/or magnetic order [|^. In the present case it is 



very difficult to determine self-consistently all the expectation values which appear. Green 
therefore approximated them by their values in the atomic limit. It then turns out that in 
the zero field paramagnetic state, for J = cxd and with the elliptic band, the CPA Green 
function G again satisfies equation ( ^.121) , with G^^ now given by equation ( |6^ . 



The densities of states calculated for T = using equations ( |6.9| ) and ( |4.12| ) with 
5* = oo, n = 0.5, u/W = 0.05 and various values of g/W are shown in figure |16|. Apart 
from lacking the perfect symmetry about the chemical potential /i = the results are 
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qualitatively similar for other values of n not too close to or 1. For g = we recover 
the elliptic band with half-width as for the DE model with J = oo, S* = oo. 

As g increases the density of states broadens and small subbands are split off from the 
band edges. As g increases further a pseudogap develops near the chemical potential. At 
a critical value g = gc a gap appears which contains a small polaron band around the 
chemical potential. Increasing g further causes more bands to be formed in the gap, with 
weights similar to those of the relevant atomic limit. It should be pointed out that the 
paramagnetic state considered here at T = is not the actual ground state, which is 
ferromagnetic. We discuss the magnetic state later. The effect of increasing temperature 
on the density of states in the gap region is shown in figure |l^ for = O.lSiy > gc- With 
increasing T the polaron bands grow rapidly and eventually merge to fill the gap. 



g/w=o 
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Figure 16: The one-electron density of states (DOS) for the Holstein-DE model with half- 
bandwidth W, for the hypothetical paramagnetic state at T = 0, with various strengths of 
electron-phonon coupling g/W . Other parameters as in figure 

mi- 



It is important to compare these results with the standard small polaron theory de- 
veloped by Holstein [p.04 |105|| . Holstein distinguished between 'diagonal transitions', in 



which the number of phonons is unchanged as the electron moves from site to site, and 
'nondiagonal transitions' in which phonon occupation numbers change. The former give 
rise to a coherent Bloch-like polaron band of half-width We~'^^'^''^^^ which decreases with 
increasing temperature. The nondiagonal transitions are inelastic processes which destroy 
phase coherence and the polaron moves by diffusive hopping. The hopping probability 
increases with temperature so that polaron motion crosses over from coherent Bloch-like 
at T = to diffusive hopping as k-^T approaches the phonon energy uj. The paramagnetic 
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Figure 17: Evolution with temperature [3 
pseudogap around the chemical potential /z 
can just be seen in figure IT6| p2|. All parameters as in figure f5. 



= {k-BT) ^ of the polaron subbands in the 
for g/W = 0.18. These subbands at T = 



state of the Holstein-DE model differs from this standard picture in one important respect. 
There are no well-defined Bloch states, owing to strong scattering by the disordered local 
spins, so no coherent polaron band will form. This is fortunate because the CPA treat- 
ment of electron-phonon scattering will never lead to coherent states of infinite lifetime 
at the Fermi surface at T = 0. However in the presence of strong spin disorder it should 



be satisfactory. We interpret the central band around the chemical potential in figure |T7 
as an incoherent polaron band whose increasing width as the temperature rises is due to 
life-time broadening of the atomic level. The life-time decreases as the hopping probability 
increases with rising temperature. 

To substantiate this picture we study the central polaron band in the limit of very 
strong electron-phonon coupling. In this limit it can be shown that we need retain only 
the r = term in equation ( |6.9D and it is then easy to solve equation ( [1.12| ) for G. The 
result is of the same form as equation (|4.15| ) but with 



= (2X [b{b + 1)]^/') (6.10) 

and = D^/W^ in the case S = oo. The central band is thus elliptical with half-width 
D and weight D'^/W'^. It is now easy to calculate the conductivity a from equation ( [4.18| ) 
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and, using D ^ W , we find 

The last step follows by using the asymptotic forms for strong coupling and high tem- 
perature Io(-2;) ~ (2^) ^^^expz and b ~ {l3uj)~^ . This form of a is similar to that for 
small polaron hopping conduction in the adiabatic limit {W ^ u) but with activa- 
tion energy Xuo/A equal to one quarter, instead of one half, of the polaron binding energy. 
Nevertheless this establishes the link between the many-body CPA and standard small 
polaron theory in the strong coupling limit. However the results shown in figure |l^, with 
parameters relevant to typical manganites, are far from this limit. They correspond to 
intermediate coupling and in the actual paramagnetic state above the Curie temperature 
the polaron bands are largely washed out. In this regime, with increasing temperature, 
there is a crossover from polaronic behaviour to a situation where the phonons behave 
classically, the case considered by Millis et. al. ||88|. For electron-phonon coupling greater 



than a critical value these authors find a gap in the density of states which gradually fills 
with increasing temperature. However in their classical treatment there are no polaron 
bands in the gap so that the link with standard polaron physics is not established. 

Apart from the symmetry about e = the above results for n = 0.5 are not untypical 
of the general case. For general n the main lower and upper bands, separated by a gap 
for g > gc, have approximate weights n and 1 — n respectively. The chemical potential at 
T = is always confined to the polaron band arising from the r = term of equation (|6.9| ), 
and moves from the bottom at n = to the top at n = 1, so that we correctly have an 
insulator in these limits. 

To calculate the Curie temperature Tq we need the full CPA theory combined with an 
exact result of DMFT for 5 = oo p2l. Results on Tc for the same parameters as before 



are plotted as functions of electron-phonon coupling g in figure [1^. The suppression of 



Tc with increasing g was first noted by Roder et. al. |^] and the CPA results are quite 



similar to those of Millis et. al. In CPA there is no reliable means of calculating the 
probability distribution function P{S^), so to go below Tc Green ||4^ used the mean-field 
approximation for the ferromagnetic Heisenberg model with classical spins and nearest 
neighbour exchange. The exchange constant is determined by Tc. We plot the up- and 



down-spin density of states for T = 0.005W/kB -C Tc and g = 0.1614^ > gc in figure |19|, 
also showing curves for the saturated ferromagnetic state and paramagnetic state at T = 
for comparison. The value g = 0.161^ is closer to gc than the value of 0.181^ used in 



figures |1^ and [l^ and we discuss these results in relation to the manganites in section In 
figure ^ the resistivity p is plotted as a function of temperature, for the same parameter 
set, with different applied fields h. The resistivity peaks sharply at Tc, and for comparison 
we show results for weaker electron-phonon coupling g/W = 0.10 in figure p^. The curve in 



figure 21 is almost indistinguishable from that of figure 7 in reference |42| for g/W = 0.01. 
This is not surprising since we see from figures |16| and |18| of this paper that the density of 
states and Tc change very little between g/W = and g/W = 0.1. These results are all 
discussed further in section ||. 
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Figure 18: Suppression of the Curie temperature of the Holstein-DE model with increasing 
electron-phonon coupling g/W. The plot is for S" = J = oo, /i = 0, n = 0.5 and u/W = 

0.05 m. 



7 Polarons and bipolarons 

Much experimental data on the manganites in both the paramagnetic and ferromagnetic 
state is interpreted in terms of the standard Holstein small-polaron theory ||104| , |105| , p9| . 
However there are usually difficulties in finding parameters which are both plausible and 
capable of explaining data from more than one type of measurement. This is not surprising 
in view of the conclusions of the last section, namely that many manganites are only just at 
the threshold of small-polaron formation even in the paramagnetic state. This means that 
one is still far from the strong coupling limit where standard theory applies. Let us review 
some of the problems in applying the standard theory, starting with the low temperature 
state. 

The width of the coherent polaron band at T = is reduced from the bare band-width 
by a factor e~^p/'^, where Ep is the polaron binding energy {g'^/to in the notation of the last 
section) and u is the phonon energy. For the typical values Ep = 0.5 eV, u = 0.05 eV this 
factor is less than lO"'*. The polaron band is so narrow that any disorder in the system will 
produce Anderson localization and a metallic state is impossible. To avoid this problem, 
by reducing the narrowing factor to 10~^ say, one requires Ep/u ~ 2 so that Ep ^ 0.1 eV. 
But then the condition for small-polaron formation Ep > Wz~^ , where W is half the bare 
band- width and z is the number of nearest neighbours ||100|| , is far from being met. Even a 
band-narrowing factor of 0.1 is too much to be compatible with heat capacity measurements 
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Figure 19: The up- and down-spin density of states of the Holstein-DE model with g/W = 
0.16 for ksT = O.OObW < ^bTc where {S^) = 0.915. Also shown are the DOS for the 
saturated ferromagnetic state at T = and for the hypothetical paramagnetic state at 



T = 0. All plots are for ^ = J = cx), /i = 0, n = 0.5, uj/W = 0.05 and g/W = 0.16 ||2 
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Figure 20: Resistivity p versus temperature for the Holstein-DE model with S = J = oo, 
n = 0.5, u/W = 0.05, intermediate coupling g/W = 0.16 and various applied fields h. 
The lattice constant is taken as a = 5 A, slightly larger than the Mn-Mn spacing in the 



manganites ||42| . 




0.00 



0.02 0.04 0.06 0.08 
Temperature k^T/W 



0.10 



0.12 



Figure 21: The same plot as figure but for weak electron-phonon coupling g/W = 0.10. 
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the term linear in T at low temperatures indicates a quasi-particle density of states 
not much larger than that given by band calculations. 

To examine the high temperature state one must consider the adiabatic limit (hopping 
integral > phonon energy). Worledge et. al. ||112|| used the result of Emin and Holstein 
||111|| for the conductivity in this limit, a = (A/T) exp{—Eg,/kBT), to obtain excellent fits to 
their data on Lai-xCaxMnOs films for T = 300-1200 K over the whole range of x. However 
for X = 0.3 the activation energy Eg, ^ 0.08 eV which is the same as in Ndo.rSro.sMnOs 
films ||113|| ; according to the theory this should be about Ep/2 which makes Ep too small 
for consistency with small-polaron formation. Also according to small-polaron theory |^ 
the optical conductivity peaks at phonon energy 2Ep ^ '^E^,] almost universally this peak 
in the paramagnetic state of pseudo-cubic manganites is observed at about 1 eV fS^ which 
is therefore inconsistent with E^. ~ 0.08 eV. This problem does not occur in the many-body 
CPA calculations where Eg_ and the optical conductivity peak are both given correctly for 
one set of parameters |^ (see section |8.3.2| ). 

Alexandrov and Bratkovsky (AB) [|114 1115| , |116|| have recently proposed a new polaronic 
theory of the manganites, including CMR. It is based on an extended Holstein model ||117 



in an attempt to avoid some of the difficulties discussed above. AB assume the O p-hole 
model of the manganite electronic structure so that DE does not come into play. The 
antiferromagnetic interaction between the local Mn spins S = 2 and the carriers is then 
assumed weak enough to treat within mean field theory. It is supposed that small polarons 
are the carriers in the ferromagnetic state, but near Tq they mostly combine to form 
immobile singlet, and possibly triplet, bipolarons. This leads to the rapid rise in resistivity 
near Tq; it falls slowly above Tq as the bipolarons dissociate thermally. It is necessary 
to assume that it is not favourable for polarons in the ferromagnetic state to combine 
as triplet bipolarons. AB propose that this is because the immobility of the bipolarons 
suppresses their exchange interaction with the local spins. The two holes in the bipolarons 
are supposed to reside on nearest neighbour oxygen sites. An applied magnetic field tends 
to split up the immobile bipolarons and reduce the resistivity; this is the CMR effect. 

To avoid the problem of very narrow polaron bands at T = Alexandrov and Ko- 
rnilovich 117 | introduced an extended Holstein model (EHM) in which the interaction 
between an electron on site i and a local phonon at site j is long range, corresponding to 
the effect of unscreened Coulomb interactions. This may also be regarded as an extension 
of the Frohlich dielectric continuum model ||118|] to a discrete lattice. The model has also 
been studied by Fehske et. al. ||119|| . They stress that the EHM polaron is a large polaron 
in the whole electron-phonon coupling range. That is the lattice distortion is spread over 
large distances even if the coupling is strong, a regime where a small polaron is formed 
in the Holstein model. Alexandrov and Kornilovich describe a 'small' Frohlich polaron, 
with a large size lattice distortion but an electronic wave function of small radius. It is 
unclear how these two length scales can be distinguished ||119|| unless the electron is tightly 
bound to a defect. The main point is that the EHM polaron band is much less narrow 
than the small Holstein polaron band with the same binding energy. This is because the 
lattice distortion undergoes smaller relative changes as the electron moves from site to site. 
AB ||116|| also find that the peak energy in optical conductivity is 27_E'p, with 7 ~ 0.2-0.4, 
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and i?a = The low observed activation energy discussed above can therefore be 

accommodated but, in a theory involving single polarons only, the inconsistency with the 
energy of the peak in optical conductivity above Tq would persist. However, according to 
AB, this peak above Tc is associated with splitting up a bipolaron, which requires a larger 
energy than 2'yEp. The theory is therefore consistent with the observed shift of optical 
conductivity towards higher energy on going from the ferromagnetic to the paramagnetic 
state. This is discussed again in section ^.3.2 . 



AB's theory presents certain problems. First, the model based on O p-hole conduction 
is open to considerable doubt, as discussed in section Secondly, it is unclear why, in a 
X = 0.3 metallic manganite with residual resistivity less than 0.1 mf2cm , the long-range 
Frohlich interaction is not screened. Chakraverty et. al. ||120| , |121|| argue that it will be 



reduced to a local Holstein interaction. Also the small enhancement of the quasi-particle 
density of states at the Fermi level, as deduced from the specific heat, is hardly compatible 
with a polaronic band at T = 0. Zhao et. al. ||122|| argue for small-polaron metallic 
conduction in the ferromagnetic state of LCMO on the basis of fitting their measured 
resistivity p{T) below 100 K to the theoretical predictions. From other data they infer 
mass enhancements of 9 and 35 in Lao.rCao.sMnOs and Ndo.7Sro.3Mn03 which, although 
small for small polarons, seem too large to be compatible with specific heat measurements. 
Another consideration is that AB's 'small' Frohlich polaron is really a large polaron. If 
a bipolaron can form it will be an extended object in which the two carriers are only 
confined within the large region of lattice distortion. Nagaev argues strongly against AB's 
theory and we refer again to the public correspondence ||123| , |116|] between him and AB in 
section |8TT| . 

Roder et. al. adopted a different method of avoiding the extreme polaronic band- 
narrowing effect. They used a variational Lang-Firsov transformation which contains 
parameters 7 and A whose role is to quench part of the dynamical polaron effect as a 
static distortion. Thus the operator s, introduced just before equation ( |6^ , becomes 
— {g/u){'yn + A) (6^ — b) and the transformation is applied to every site in the lattice. The 
band-narrowing factor is then modified from e"'^'"^^ to e""^ ^^^^ and < 7 < 1. The main 
result is that the Tc of the pure DE model is reduced by this factor due to electron-phonon 
coupling; the resistivity was not calculated. The most recent work along similar lines is 



that of Perroni et. al. [124 |125|| . The DE effect is treated by a simple band-narrowing but 
these authors go much further than Roder et. al. In particular they calculate a one-electron 
Green function which is quite similar to Green's |^] and, like his, yields the correct atomic 
limit (equation (|6.9| )) for the classical spin case considered. However instead of a smooth 
evolution of a pseudogap with increasing temperature, they find phase separation above 
0.65Tc with ferromagnetic regions of low electron density (large polaron behaviour) and 
paramagnetic insulating regions of high electron density (localized small polarons). How- 
ever the effect of the nanoscale structure of such regions, enforced by Coulomb interaction, 
is not considered. 
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8 Theory and experiment 



The experimental evidence and theoretical arguments reviewed in previous sections point 
to the Holstein-DE model as a plausible model of the manganites Ai„xA^Mn03, at least 
for X ~ 0.33. For this degree of doping with the divalent element A' one is normally in the 
middle of the ferromagnetic regime and inhomogeneity due to charge and orbital ordering 
is least likely. In order to confront theory with experiment effectively it is essential to keep 
the number of adjustable parameters to a minimum. Ideally one would like to understand 
the difference between materials, as reflected in the results of many different types of 
experiments, by changes in one crucial parameter. The many-body CPA approach to the 
Holstein-DE model gives us the opportunity to attempt this and from section ^ it is clear 
that the crucial parameter is the electron-phonon coupling g/W. In this section we take 
flxed reasonable values of the parameters u and W. Band calculations [|1^, |T8| suggest 
= 1 eV as an appropriate half-bandwidth for the eg band. Then it is reasonable to 
take uj/W = 0.05 to correspond to the observed transverse optic phonons with u ~ 40-70 
meV which couple strongly to the electrons in LCMO ||126|| . It is also reasonable to take 



the DE limit J ^ oo ||127|] . Furthermore we see from figures |14(a)| and that for the 



DE model neither Tq nor the resistivity p vary enormously with S so that = cx) is a 
convenient approximation to the 5 = 3/2 Mn spins. These are the parameters which have 
been used in figures |1^ The band-filling is also restricted to n = 1 — x = 0.5 rather 
than n = 0.6-0.7. This is very convenient because the chemical potential remains fixed for 
all T by electron-hole symmetry in the case S = oo. This will change the critical value of 
g/W for small polaron formation somewhat but the correct general picture should emerge. 
Since we consider a homogeneous state we are not concerned with the existence of charge 
ordering for x = 0.5. We now compare the theory with different types of experiment. 



8.1 Transport properties and Curie temperature 

Perhaps the most striking feature of the manganites is the very different behaviour observed 
in apparently similar materials such as LSMO and LCMO. For LSMO, with x ~ 0.33, 
Tq ~ 370 K whereas for LCMO, with a similar x, Tq ^ 240 K. The difference in behaviour 
of the resistivity p above Tq is much more striking (see figures ^ and For LSMO 
p ^ 4mf2cm and increases slowly with temperature as in a poor metal ||l5l. The p{T) 



curve is very similar to that of figure |2T] for g/W = 0.1 except for a much larger resistivity 
at low temperature in the calculations. Since this feature persists even for g/W = 0.01 
(figure 7 in reference |^) it presumably arises from overestimated spin disorder scattering 
at low temperatures due to use of the classical spin Heisenberg model to determine P{S^). 
In LCMO the resistivity rises to a maximum at Tq of about 40 mf2cm and then falls with 
increasing temperature above Tq- In contrast to LSMO there is thus a transition from 
metallic to insulating behaviour. Also the resistivity peak is strongly reduced and shifted 
to higher temperature with increasing applied magnetic field (figure ^. This is the CMR 
effect. This type of behaviour is seen in figure ^ for g/W = 0.16. The decrease of p with 
increasing temperature above Tq is associated with the increasing density of states at the 
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Fermi level, as seen in figure [T^. The main differences between theory and experiment 
are a more rapid observed drop in p with decreasing temperature below Tq and a more 
sensitive observed CMR effect. h/W = 0.004 corresponds to a field of about 20 T for 
W = 1 eY and the corresponding reduction in p in figure is achieved with a field of 
about 5 T experimentally. Millis et. al. ||8^ noted a similar problem in their work using 
classical phonons. Both of the discrepancies mentioned might be remedied by introducing 
a dependence of g on p, corresponding to more efficient screening of the electron-phonon 
interaction with increasing density of states. The huge reduction in resistivity peak on 
reducing g/W from 0.16 to 0.10 shows the extreme sensitivity of p to changes in g. The 
main point to notice is that we can understand the enormous difference between LCMO and 
LSMO within the Holstein-DE model by assuming the electron-phonon coupling changes 
from g/W = 0.16 in LCMO to g/W = 0.10, or slightly greater, in LSMO. The observed 



ratio of the Curie temperatures, slightly less than 2, is then in accord with figure 18 



As discussed in section || the critical coupling gc for the formation of a polaron band is 
gc/W ~ 0.15, with phonon energy u/W = 0.05, and to obtain the right order of magnitude 
for p above Tq in LCMO, g/W is pinned down closely to 0.16. A larger value for g leads to 
too high a resistivity and too low a Curie temperature. It is interesting that neither p(T) 
nor Tc change when g/W is varied between 0.1 and 0. This means that the resistivity of 
LSMO can be described quite well by the pure DE model, as stressed by Furukawa ||6T| , 
but electron-phonon coupling is not negligible and shows up in the optical conductivity, for 
example, which we discuss in section |8.3.2| . However, from the results of Millis et. al. |Q 
for classical phonons, one can understand why a coupling small enough to give a LSMO- 
like p(T) curve does not lead to a change in slope of the rms oxygen displacements, as 
a function of temperature, at Tc- No such change is found in LSMO [H, in contrast to 
the case of LCMO ||128|| . It is more difficult to understand the observation ||129|| of static 
local Jahn- Teller distortions in LSMO at room temperature, apparently associated with 
localized carriers in the presence of metallic conduction. 

From figure |16| we see that for g/W = 0.1, appropriate to LSMO, there is no sign of a 
pseudogap in the density of states. An actual gap in the hypothetical paramagnetic state 
at T = appears at g = gc with gc/W between 0.15 and 0.16. From figure |19| we see that 
for g/W = 0.16, appropriate to LCMO, a few polaron subbands have appeared in the gap. 
These are seen much more clearly in figure for g/W = 0.18 when there is a larger gap. 
However the subband structure is washed out completely for PW = 25, corresponding to 
T = 464 K for ly = 1 eV, and this effect will occur at a much lower temperature for 
g/W = 0.16. Thus in the actual paramagnetic state of LCMO above Tc we do not expect 
the quantum nature of phonons to manifest itself, so we are essentially in the classical 
regime of Millis et. al. |^]. The same is true in the saturated ferromagnetic state at 
T = where only a pseudogap appears in figure |l^. As the temperature rises towards 
Tc a minority spin band grows, also with a pseudogap, while the majority spin band loses 
weight. The width of the ferromagnetic bands decreases with increasing temperature, 
corresponding to the DE effect, but the narrowing in the paramagnetic state is not so 
marked as in the pure DE model. Thus double exchange is not so effective in the presence 
of strong electron-phonon coupling, which is consistent with the reduction in Tc shown in 



49 



figure 0. We discuss this again in section ^]4| in connection with spin-waves. 

Clearly the picture of the manganites which emerges here is close in spirit to that of 
Millis et. al. ||8^, although the relationship to polaron physics is not so clear in their clas- 
sical approximation. We have described earlier how some other authors adopt completely 
different viewpoints. In particular Nagaev ||130|| argues against any polaronic effects, while 
Alexandrov and Bratkovsky ||115| , |116|| assume strong electron-phonon coupling with small 



polarons even in the ferromagnetic state and with immobile bipolarons forming near Tq. 



The public correspondence ||123|, 116| between Nagaev and Alexandrov and Bratkovsky 



(AB) centres on estimating the magnitude of the polaron binding energy Ep and the cri- 



terion for small polaron formation ||114|| . Since the picture of LCMO presented here lies 
between their extreme views it is interesting to compare our estimates with theirs. For 
LCMO we find Ep = g^/u ~ 0.5 eV ioi W = 1 eV whereas, for manganites in gen- 
eral, Nagaev estimates Ep ^ 0.1-0.3 eV and AB estimate i?p ~ 1 eV. Our condition for 
small-polaron formation in a paramagnetic state at T = is > (7c ~ 0.1514^ which corre- 
sponds to Ep > 0.45iy. Nagaev adopts the criterion Ep > W, remembering that W is the 
half-bandwidth in our notation, while AB ||114|| propose Ep > 2iy(82)-i/2 _ Q 29W with 



number of nearest neighbours z taken as 6. AB's condition is less stringent than Eagle's 
lOOll condition for 'nearly small polarons' Ep > Wz~^l'^ = OAIW which is close to ours. 



Both as regards this criterion and the value of Ep for LCMO, the results of the Holstein-DE 
model are intermediate between Nagaev's and AB's, as expected. For LSMO, on the other 
hand, our estimate of Ep is 0.2 eV. In this case we agree with Nagaev and Furukawa |^ 
that electron-phonon coupling is not so important. AB's work is reviewed in section ^ 

8.2 Isotope and pressure effects 



There is an excellent review of oxygen isotope effects in manganites by Zhao et. al. ||131 
The most striking effects are a decrease in Tq and an increase in resistivity p, partic- 
ularly near Tq, when ^^O is replaced by ^^O. However other important effects are an 
induced metal-insulator transition, a shift of the charge-ordering transition and effects on 
the thermal-expansion coefficient and electron paramagnetic resonance (EPR) measure- 
ments. Although Zhao et. al. attempt to understand these effects using the concepts of 
small-polarons and bipolarons, they point out clearly the difficulties in such a theory. We 
believe that most of these difficulties can be resolved within the many-body CPA approach 
to the Holstein-DE model, which frees us from the strong coupling regime of standard 
small-polaron theory. 

The oxygen isotope exponent is defined as ao = — (ATc/Tc)/(AM/M), where Tq and 
M are the Curie temperature and oxygen mass for the ^^O sample, ao ranges from 0.07 
with Tq = 367 K in Lao.eySro.asMnOa to 0.85 in Lao.8Cao.2Mn03+y {Tq = 206 K) and up to 



4 in (Lao.25Ndo.75)o.7Cao.3Mn03 {Tq ~ 110 K). Zhao et. al. ||131|| point out that all available 
data (for x ~ 0.2-0.35) can be fitted by an equation 

ao = 21.9exp(-0.016Tc) . (8.1) 
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According to the DE model and small-polaron theory Tc oc exp(— 7^£'p/co'), which is the 
usual band-narrowing factor (see section |^) with < 7^ < 1 [^. When the mass M 
changes, u varies as M~^/^ but the polaron binding energy is independent of M since 
it corresponds to a static response of neighbouring oxygen atoms to an eg electron on 
a Mn site. Hence Tc decreases rapidly with increasing M and it is easy to show that 
Tc oc exp(— 2a;o) ||131J] . Zhao et. al. ||131J] point out that this is incompatible with the 
experimental relation ( ^.11) . They go on to make the very important observation that the 
effect of pressure on Tc, again over a wide range of manganites, is summarized accurately 
by the equation 

dlnTc/dP = 4.4exp(-0.016Tc). (8.2) 

Combining equations (|8.1| ) and (|]2|) one has ao = 5(dlnTc/dP). Zhao et. al. deduce 
from this simple proportionality between the isotope exponent and the pressure-effect co- 
efficient that the major effect of pressure is to increase the phonon frequency. We come 
to a completely different conclusion from the following analysis using the many-body CPA 
approach to the Holstein-DE model. 



In discussing the pressure effect, Green |^ observed that the increase in Tc, and de 



crease in p, stems mainly from a decrease in the effective coupling constant g"^ / {uW) = 
Ep/W. As pointed out in section |l| (see figure |), the effect of pressure is equivalent to an 
increase in (ta) and broadens the eg band, that is W increases. Green |^ therefore mod- 
elled the strong suppression of the resistivity peak and the increase in Curie temperature 
in LCMO by increasing the bandwidth, and assuming other terms in the Hamiltonian are 



constant. Calculated results are shown in figure |22|. A simple estimate, using the known 
compressibility and dependence of W on lattice constant, shows that the theoretical pres- 
sure for a given effect is about four times larger than that required experimentally. This is 
the same factor that was found in the case of the magnetic field required for a given CMR 
effect, so both discrepancies could possibly be removed with the same dependence of g on 
p, due to screening, postulated in section |8T| . We return to this point later. In fact near Tc 
the CMR effect (figure ^) and the pressure effect (figure ^) are both driven by a change 
of bandwidth. In the first case the magnetic field produces a substantial magnetization for 
T ^ Tc and this increases W by the DE effect. 

To investigate the isotope effect in our model, assuming only an effect via u, some care 
is needed. It was pointed out that in equation ( |6.1|) the electron-phonon coupling term 
corresponds to a term —g' riiXi where Xi is to be associated with oxygen displacement 
around a Mn atom. Here g' should be independent of the oxygen mass M. In the second- 
quantized form of equation ( |6.1J ) one finds g = g' {2Muj)~^^'^ so that the polaron binding 
energy Ep = g'^/ui = g''^ / (2Mu;^). Since for an oscillator 00 oc M~^/^ the polaron binding 
energy g'^/co is independent of M as expected. However g varies as M~'^^^. We therefore 
recalculated the resistivity p(T) with the same parameters as used in figure ^for ^^O but 



with g and u scaled appropriately for ^^O. The results are compared in figure |2^. The 
almost complete absence of an isotope effect in both Tc and p confirms Green's observation 
that p(T) depends almost exclusively on the parameter Ep/W. Since Ep should not depend 
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Figure 22: The effect of pressure on the resistivity and Curie temperature in LCMO by 
increasing the half-bandwidth W to W = 1.05W and I.IW li^ . 



on M we have to conclude that the isotope effect is due to a decrease of bandwidth 2W 
when M is increased. 

Thus we make the hypothesis that the isotope effect occurs via the eg bandwidth 2W 
rather than the phonon energy u, whose undoubted change has no effect. We must stress 
that 2W is the bare band-structure bandwidth. Clearly a sort of scaling theory, linking 
CMR, pressure effect and isotope effect, can be constructed but this will be done elsewhere. 
The change in the p(T) curve is typically quite similar for ^^O ^^^O back-substitution to 
an applied field of 5-10 T or a pressure of 5-10 kbar. Babushkina et. al. ||132|| have 
already suggested that a metal-insulator transition induced by oxygen isotope substitution 
in (Lai_yPry)o.7Cao.3Mn03 may proceed by a change of bandwidth. They point out that 
the mean-square displacement of an ion from its nominal position depends on the ionic 
mass M, even at T = 0, due to zero-point vibrations. Such displacements could affect 
the Mn-O-Mn bond angle on which the bandwidth critically depends. Similar effects 
were observed earlier in (Lao.5Ndo.5)o.67Cao.33Mn03 ||133|| . Babushkina et. al. also see large 
shifts of the p(T) curves with applied magnetic field, as one would expect from the above 
discussion. Clearly similar ideas can be applied to shifts in charge-ordering transitions due 
to isotopic substitution and magnetic field. The isotope effect in the jump of the thermal 
expansion coefficient at Tq is easily understandable since it is related by thermodynamics 



to the pressure effect 131 



We now have to address the question of whether replacing ^^O by ^^O can produce the 
required decrease of W to account for the effect on p{T) . The change in p{T) is typically 
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Figure 23: The absence of an isotope effect in LCMO due to a simple mass scaling 
(oc M~^/^) of the phonon frequency. 



like that calculated for a 5-10% decrease in W (cf. figure ^). If we allow for the factor 
of 4 proposed earlier, due to improved screening of the electron-phonon coupling as the 
system becomes more metallic this is reduced to 1-2%. According to Zhao et. al. ||134|| the 
volumes of the unit cells of ^^O and ^^O samples of Lao.eTCao.ssMnOs at room temperature 
are the same within the accuracy of the x-ray determination. Singh and Pickett ||23| 
have investigated the change in band structure of a La2/3Bai/3Mn03 virtual crystal on 
going from a cubic perovskite structure to a Pnma structure consistent with the neutron 
refinement of Dai et. al. ||128|| for Lao.esCao.ssMnOa, the two structures having the same 
unit cell volume. Large band-narrowings occur owing to bending of the Mn-O-Mn bonds; 
in particular a gap of a few tenths of an eV opens between the majority spin eg and t2g 
bands. The distortion from the cubic structure has no static JT components, consisting 
essentially of rotations of O octahedra. It is not inconceivable that isotopic substitution 
could modify these rotations sufficiently to give a 1-2% reduction in the eg bandwidth. 
Further theoretical and experimental investigation of this possibility is highly desirable. 
Finally we come back to the enhancement factor of 4 postulated above. It means that 
a 1% increase in W corresponds to a 4% decrease in Ep/W, the large relative change in 
Ep being assumed to arise from increased screening with increased metallization. This is 
not unreasonable since it is known from optical conductivity measurements, discussed in 
section |8.3.2| , that the polaron binding energy Ep decreases from 1.2 eV in Ndo.ySro.sMnOs 
to about 0.6 eV in LSMO. The change in Cg bandwidth between these two materials is 
certainly much less than a factor 2 so that in general a much larger relative change in 
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Ep/W than in is a reasonable hypothesis. 

In summary our picture of the isotope effect due to ^^O — i>^^0 substitution is as follows. 
The change in phonon frequency plays no part. The effect is driven by a 1-2% reduction 
in the eg bandwidth 2W associated with small changes in the rotations of O octahedra 
probably arising from modified vibrational amplitude of the O ions due to their mass 
change. This primary effect is strongly enhanced by an increased polaron binding energy 
Ep arising from reduced metallic screening. The scaling parameter Ep/W links CMR, 
pressure-effect and the various isotope effects. 

We should mention that Nagaev [|130| , |135|| has taken a completely different line on 
the isotope effect in which electron-phonon coupling has no role. He considered several 
possibilities, all based on the isotope dependence of the number of excess or deficient 
oxygen atoms in thermodynamic equilibrium. Since he quoted experimental evidence for 
this ||136|| his arguments were quite persuasive. However Zhao et. al. ||137|| have criticized 



the experimental work of Franck et. al. [|136|| and they refute Nagaev's model in which the 
isotope effect is not intrinsic. 

8.3 Pseudogaps 

According to theory we expect pseudogaps at the Fermi level to be observable in the density 
of states of LCMO both below and above Tq- These should appear in experiments such as 
scanning tunnelling spectroscopy, photoemission and optical conductivity measurements. 
No pseudo-gap is expected in LSMO below Tc and the effect should be much less than in 
LCMO above Tc. The pseudogap is a feature of the atomic limit, typified by the envelope 
function in figure |15| for T = with maxima determined by the polaron binding energy 
g'^/w, and is completely washed out when g'^/^uW) < 0.2 (see figure [1^). Early results of 
scanning tunnelling spectroscopy on LCMO ||138|| with x = 0.3 seem unlikely to relate to 



the bulk. In the ferromagnetic state at 77 K there is a huge gap of about 1 eV. It is not clear 
why the authors interpret this as evidence for half-metallic ferromagnetism. A gap of this 
size associated with small-polaron formation in the bulk would imply an unrealistically 
large electron-phonon coupling, certainly incompatible with metallic conduction and a 



Curie temperature of reasonable magnitude for LCMO. More recently Biswas et. al. ||139 
reported a scanning tunnelling spectroscopy study of several manganites. The results are 
very much in accord with the theory. There is no gap in the low temperature ferromagnetic 
state but a small gap (pseudogap) appears for T Tc in the low (high) Tc materials. As 
T increases above Tc the pseudogap or gap gets filled in as we would expect (see figure [l^). 
The decrease in resistivity with increasing T is due to the gradual filling of the pseudogap. 

8.3.1 Angle-resolved photoemission spectroscopy (ARPES) 



In an extremely interesting paper on ARPES for the bilayer manganite Lai.2Sri 8Mn20 



7, 



nominally with n = 0.6, Dessau et. al. ||36|| interpret their results very much in the spirit of 
the Holstein model. The low Tc = 126 K in this bilayer manganite is partly due to quasi- 
two dimensional fluctuations, but the large resistivity p ~ 3 mQcm at low temperatures 
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indicates that small-polaron bands might exist even in the ferromagnetic state. Conse- 
quently the electron-phonon coupling should be stronger than in cubic manganites like 
LCMO and to model the system by the Holstein-DE model Hohenadler and Edwards |Q 
chose g/W = 0.2. The one-electron spectral function Ak{e) is given by equation (|4.17| ) and 



the band is taken to be of the form = ~W cosny for k = 7r(l,|/), <y <1 with W = 1 



eV as usual. This roughly models a band, calculated by Dessau et. al. ||3§|, which crosses 
the Fermi level i^p at = 7r(l, |). The calculated results for are shown in figure |24|. 
Well away from the Fermi level, a well-defined peak exists which broadens as k approaches 
the Fermi momentum at y = 0.5. For larger y the weight below the Fermi level is strongly 
reduced. The peaks never approach the Fermi level closely which is an important feature 



of the observed spectra [^] reproduced in figure |25|. The theoretical curves in figure |24 
resemble quite closely the data of figure |25|(c). There is a pseudogap in the calculated 
spectra extending about 0.1 eV on each side of the Fermi level. In fact this pseudogap 
contains polaron bands like those shown in figure O. However, their amplitude is too 



small to show up in figure ^ and in the experimental data. Nevertheless, it is the central 
polaron band around the Fermi level which is presumably responsible for the low but finite 
conductivity of the system. This comparison between theory and experiment supports the 
conclusion of Dessau et. al. |^ that, in the manganites with a layered structure, strong 
electron-phonon coupling (with the appearance of a pseudogap) is already important be- 
low Tc- This contrasts with the usual pseudocubic manganites where the pseudogap only 
appears above Tc. It should be mentioned that Moreo et. al. ||140|| interpret the observed 



pseudogap not as an intrinsic property but in terms of phase separation. 
8.3.2 Optical conductivity 



Using the notation of section we can write the optical conductivity cr(z/), corresponding 
to the elliptic density of states in equation ([4.8|), in the form ||141| , |142| , ^ 

cr(//) = dej dEM,{E)AE{e)AE{e + y)M^JS^±^ . (8.3) 



This expression satisfies the correct one-band sum rule ||142|, |50| that (2/7r) a{u)du 



-Ke^/{3ah), where the 'kinetic energy' K is the thermal average per lattice site of the 
first term in the Hamiltonian ( |6.1|) . The energy K appeared previously in equation ( p.ll|) 
for the spin- wave stiffness constant. For = 0, equation (|8.3|) yields the dc conductivity 



given by equation ( |1.19|) . In the ferromagnetic state of the pure DE model {g/W = 0) at 
T = the spectral function AE^e) = 6{e — E) so that, from equation (|8.3|) , a{u) oc 5{u). 
However, in LSMO with x = 0.3 the observed intraband transitions spread up to 1 eV ||143| 



which suggests that phonons are important for producing broadened spectral functions. A 
sharp Drude peak below 0.04 eV, arising from 5-like parts of AE^e) for E near the Fermi 
level, accounts only for half of the total intraband spectral weight [|143|| . An alternative 
interpretation is that orbital degrees of freedom in the doubly-degenerate eg band, combined 
with strong correlation (no doubly-occupied sites), lead to incoherent motion of carriers 
144], |145||. Here we pursue the comparison of results of the many-body CPA treatment of 
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Figure 24: The spectral function AyJ^e) in the ferromagnetic state at T = for the Holstein- 
DE model with J = S* = oo, n = 0.5 and strong electron-phonon coupling gjW = 0.20, 
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T = lOK, reproduced from Dessau et. al. p^j . 
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the Holstein-DE model with experiment. A defect of Green's |^ CPA treatment is that the 
sharp quasi-particle peak in the spectral function, which should exist in the ferromagnetic 
state at T = 0, is missing. This is due to incoherent scattering by the phonons which gives 
a spurious residual resistivity. Thus the Drude peak in cr(z/) is absent, as in the classical 
phonon treatment of Millis et. al. [^]. Also in the one-band model with J = oo interband 
transitions between Hund's rule split bands, and p— >d charge transfer transitions do not 
feature. These are observed at photon energy u > 3 eV and cr(z/) calculated in the one-band 
model is only non-zero in the region u < 2.5 eV with in general only one peak. Hohenadler 
and Edwards calculated cr(z/) in the Holstein-DE model with g/W = 0.16 and compared 
with experimental data on an unannealed film of Ndo.ySro.sMnOa (NSMO) ||146|| . The 



comparison is shown in figure |26|. This data was chosen because the low temperature 
incoherent scattering in the sample matches closely the incoherent scattering introduced by 
CPA. In annealed NSMO films [^], and in single crystals [|1471| , cr(z/) in the low temperature 
ferromagnetic state continues to rise with decreasing u down to much lower photon energy, 
and o"(0) ~ 3 (mficm)"^ . Clearly in our calculation, and those of Millis et. al 



the shift of the peak to lower energy is held up due to spurious incoherent scattering in 
the ground state, which limits the low-temperature dc conductivity. There seems little 
doubt that a more correct treatment of the Holstein-DE model at low temperatures would 
lead to something more like the shift observed in good samples. Very recently Perroni 
et. al. [|124|| seem to have achieved this together with a Drude peak at low temperatures in 
their variational Lang-Firsov treatment. It is to be hoped that these features survive even if 
their predicted phase separation is suppressed by Coulomb interaction. In any case one does 
not need to invoke a change from unbinding bipolarons to unbinding polarons to explain 
the shift, as is done by Alexandrov and Bratkovsky ||114|| . It is also not clear why in the 
bipolaron regime these authors find a Gaussian form for a^u) without a threshold photon 
energy. In the bipolaron regime there should be a gap in the single-particle spectrum, 
analogous to that in a superconductor, but with magnitude of the order of the bipolaron 
binding energy. There is no such gap in the metallic 'unpaired' single polaron regime. 

In the paramagnetic state above Tq, <j{i^) is much less sample dependent and the 
CPA calculations are much more reliable. Many authors ||88|, |5^, |146| , ^] agree on the 
following interpretation of the peak at about 1 eV. During the absorption process an 
electron moves from one site to a neighbouring one which was previously unoccupied. 
The electron motion is accompanied by a lattice distortion, of Jahn-Teller type, which 
corresponds to a displacement of the local phonon oscillator coordinate in the Holstein-DE 
model. When an electron enters (leaves) a site the final displaced (undisplaced) oscillator 
is generally in an excited state with typical excitation energy g'^/u. This is the atomic- 
limit polaron binding energy and for the parameters assumed here is about 0.5 eV. Thus 
the peak in a^u) (figure |26(a)|) occurs at about twice the polaron binding energy just as 
in the standard small-polaron theory ||9^. However for the present intermediate electron- 
phonon coupling g/W = 0.16 polaron bands near the Fermi level are largely washed out 
above Tq [0, so standard small-polaron theory is not expected to apply to cr(z/) for low 
photon energies. This is certainly the case because small-polaron theory would predict an 
activation energy in the dc conductivity a{0) of 0.25 eV, half the polaron binding energy. 
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Figure 26: (a) Calculated optical conductivity for strong electron-phonon coupling g/W = 
0.16 in the ferromagnetic state at T = (full line), the paramagnetic state at T = Tc 
(dotted line) and the paramagnetic state at T = 1.5Tc (dashed line). The plot is for 
J = S* = oo, n = 0.5 and a = 5 A ||4^. (b) optical conductivity of Ndo.ySro.sMnOs at 
different temperatures, reproduced from Kaplan et. al. ||146|| . 



However in NSMO the activation energy is observed to be about 0.08 eV ||113|| and Green's 
[i^ calculations (see figure ^Up are in good agreement with this. 

We conclude that the energy of the peak in a^u) for T > Tc is a reliable measure of 
twice the polaron binding energy, but that the activation energy for dc conductivity gives 
no such direct information. The most extreme illustration of this is LSMO with x = 0.3. 
There is no activation energy for cr(0) = which slowly decreases with temperature 
above Tc (see figure H). But the temperature-dependent part of a{h') has a clear peak at 
about 0.6 eV for T > Tq ||143|| . This indicates a polaron binding energy of 0.3 eV, which 
corresponds to g/W ~ 0.12 for the other parameters assumed in this section. The estimate 
of ~ 0.1 we have sometimes suggested for this material is probably too low. However 
the difference between g/W = 0.16 for LCMO and NSMO and g/W = 0.12 for LSMO, all 
with X = 0.3, is sufficient to explain their very different behaviour. 



8.4 Spin waves 



In section ^]2| we pointed out that the simple result of equation ( |3.11| ) for the spin-wave 
stiffness constant D applies even when local electron-phonon coupling is included. In 
particular it applies to the Holstein-DE model and a simple variational argument shows 
that is is an upper bound to the true value of D [^ . The complete spin- wave dispersion 



curve takes the form of that for a simple nearest-neighbour Heisenberg model p4|, H3|, as 
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observed in Lao.TPbo.sMnOs [|148|| . Deviations from this simple result are found near the 



zone boundary in manganites with lower Tq than the Pb doped material ||149| , |150| , |151| 



Hohenadler and Edwards |^ calculated D as a function of the electron-phonon coupling 
g/W and the result is shown in figure [27|. This behaviour of D in the Holstein-DE model is 
very similar to that of Tc, as calculated by Green |Q (see figure [TSl). The main difference is 
in the extreme strong-coupling limit where Tc becomes very small at g/W ~ 0.35 whereas 
D is decreasing quite slowly. The slow decrease of D is exactly what one expects from 
equation and small-polaron theory, where the kinetic energy K ~ g'"^ ||152| , 



153 . T, 



seems to be determined more by the width of the narrow polaron band around the Fermi 
level, which decreases exponentially with g. This is the result found by Roder et. al. |^ 
which is essentially Tc oc D{g = 0) exp{—'yEp/u), in the notation of section |^. Here 
D{g = 0) is given by equation (|3.11j ) for the pure DE model. For manganites with higher 
Tc, where the spin- wave dispersion curve is Heisenberg-like, one expects that the ratio 
S = D/{k-BTca?), where a is the lattice constant, should be approximately 0.286 which 
is the value for the S' = 3/2 Heisenberg model. In fact Hohenadler and Edwards |Q 
find, combining their results with those of Green for Tc (see figure |18D, that 5 ~ 0.24 for 
g/W = 0.1 which is appropriate to a high Tc manganite like LSMO {x ~ 0.3). Clearly, from 



figures |18| and ^ 6 increases as g/W becomes larger, that is as Tc becomes smaller. This 
agrees with experiment; for example in Ndo.rSro.sMnOs (Tc = 197.9 K), 5 ~ 0.64 [|149|| . 
In this system, and in LCMO ||154|| , the spin- wave stiffness constant does not collapse 
to zero at T = Tc so the behaviour is certainly not Heisenberg-like. The properties of 
the Holstein-DE model in this parameter regime need further investigation. The phase 
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separation for T > 0.65Tc predicted by Perroni et. al. ||124|| should be examined in the 
presence of Coulomb interaction. 



9 Conclusions 

A brief review inevitably involves a personal choice of topic and emphasis, Here we have 
concentrated on the physics of manganites such as Lai-xCaxMnOs in the ferromagnetic 
regime with x ~ 0.2-0.4. We have not discussed interesting phenomena such as charge 
and orbital ordering which occur elsewhere in the phase diagram. 

The main conclusion is that the Holstein-DE model is capable of describing semi- 
quantitatively a wide range of experimental data on these materials. The model combines 
Zener's double exchange mechanism for ferromagnetism with the possibility of polaronic 
effects associated with strong electron-phonon coupling. These are the essential ingredients 



proposed by Millis et. al. ^]. The double-exchange mechanism relies on the mobile 
carriers in doped systems being Mn eg electrons rather than O p holes. The theoretical and 
experimental evidence that this is indeed the case is reviewed in section ||. In section || the 
concept of double exchange is set within the general scheme of itinerant electron magnetism. 
One ingredient omitted in the Holstein-DE model is A-site disorder and the discussion of 



section |5T| shows that this omission is not serious for systems with x ~ 0.3. 

The most complete theory of the Holstein-DE model, capable of dealing with quantum 
spins and phonons, is that of Green it is developed in a concise way in section |[ The 



theory makes use of the many-body CPA method introduced earlier W^, for treating 



the pure double-exchange model, as described in section ^. For the pure DE model this 
method is equivalent to Furukawa's dynamical mean-field theory in his limit of classical 
spins. Theoretical results for the Holstein-DE model are compared with experiment in 
section ^ One striking conclusion is that the manganites fall in the regime of intermediate 
electron-phonon coupling, just on the verge of small-polaron formation. Thus in general the 
standard results of small-polaron theory do not apply. However even in the intermediate- 
coupling regime the atomic-limit polaron binding energy Ep is still a significant parameter 
and a peak in the optical conductivity a^u) occurs at photon energy 2Ep, just as in standard 
small-polaron theory. But the activation energy for electrical conduction above Tc turns 
out to be much smaller than that predicted by small-polaron theory, in agreement with 
experiment. The key parameter in the theory is Ep/W , the ratio of the polaron binding 
energy to the half-width of the eg band. This determines the Curie temperature Tc and 
the form of the resistivity p(T); there is essentially no explicit dependence on the phonon 
frequency u in the intermediate-coupling regime. Thus the oxygen isotope effect in Tc and 
p(T) is not a consequence of the change in u associated with a change of isotopic mass. 
The observed proportionality between the isotope exponent and the pressure-coefficient 
indicates that the shift of Tc arises from a change of Ep/W in both cases. The primary 
effect may be regarded as a change in W due to a change in Mn-O-Mn bond angle, 
and this effect is strongly enhanced by a change in Ep due to modified screening of the 
electron-phonon coupling arising from a changed density of states at the Fermi level. It 
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is proposed, following Babushkina et. al. ||132||, that a change in Mn-O-Mn bond angle 



arising from isotopic substitution is due to a modified vibrational amplitude of the O 
ions owing to the mass change. This proposal, and that concerning the screening of the 
electron-phonon coupling, requires further theoretical and experimental investigation. The 
large isotope effect, pressure effect, colossal magnetoresistance and the striking difference 
between Lao.TSio.sMnOa and Lao.rCao.aMnOs are all due to extreme sensitivity to the 
parameter Ep/W in the intermediate coupling regime {E^/W ~ 0.5). 

Another point which requires investigation is the nature of the ferromagnetic transition 
at Tq. It is found, both theoretically and experimentally, that in manganites with low Tq 
the ratio of spin-wave stiffness constant D to Tc is much larger than the Heisenberg-like 
value occurring in those with high Tc. Experimentally, in the low Tc systems, D does 
not collapse to zero at Tc. Nanoscale phase separation near and above Tc is a possible 
explanation [|154|| and further theoretical work along these lines is required. 

The physics of the manganites involves a subtle interplay between magnetism and 
electron-phonon coupling. The same might be said, with less certainty, of high temperature 
superconductivity in the cuprates. It seems clear from this review that a firm theoretical 
understanding of the manganites is within reach. Methods of handling electron-phonon 
coupling in the intermediate coupling regime have been developed. An intriguing prospect 
is that this experience gained with the manganites may lead to important advances in our 
understanding of the cuprates. 
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